1 64000-1 9-F 


NASA 


E86-10025 

eg- nbAbf 


(E86-10025 NASA -CB- 176269) STODY OF 
SPECTRAL/RADICMETBIC CHARACTERISTICS OF IHE 
THEMATIC MAPPER FCR LAND USE APPLICATIONS 
Final Report, 21 Sep. 1982 - 20 Jun. 1985 
(Environmental Research Inst, ci Micl^^^O 


N86- 1 6690 
THRO 

N86- 16694 
Unclas 

G3/43 00025 


Final Report 


STUDY ON SPECTRAL/RADIOMETRIC 
CHARACTERISTICS OF THE THEMATIC 
MAPPER FOR LAND USE APPLICATIONS 


21 September 1 982 — 20 June 1 985 


WILLIAM A. MALILA 
MICHAEL D. METZLER 


SEPTEMBER 1985 


C-r’ -Inal photography may ho purchat 
-o EP.CS Data Center 
3ioajc Fails , SD 

Original' photography may bo purchase 
from EROS Data Center 
Sioux Falls,. SQ 57193 

Contract NAS5-27346 
NASA Goddard Space Flight Center 
Green belt Road 
Greenbelt, Maryland 20771 


: 


RE 


CTi 


|V\$J 

40 $ 


ENV IRONMENTAL 


RESEARCH INSTITUTE OF MICHIGAN 


BOX 8618 • ANN ARBOR • MICHIGAN 48107 


NOTICES 


S ponsorshi p . The work reported herein was conducted by the 
Environmental Research Institute of Michigan under Contract NAS5-27346 
for the National Aeronautics and Space Administration, Goddard Space 
Flight Center, Greenbelt, MD 20771. Contracts and grants to the 
Institute for the support of sponsored research are administered 
through the Office of Contracts Administration. 

Disclaimers . This report was prepared as an account of 
Government sponsored work. Neither the United States, nor the National 
Aeronautics and Space Administration (NASA), nor any person acting on 
behalf of NASA: 

(A) Makes any warranty expressed or implied, with respect to the 
accuracy, completeness, or usefulness of the information, 
apparatus, method, or process disclosed in this report may 
not infringe privately owned rights; or 

(B) Assumes any liabilities with respect to the use of, or for 
damages resulting from the use of any information, 
apparatus, method, or process disclosed in this report. 

As used above, "person acting on behalf of NASA" includes any employee 
or contractor of NASA, or employee of such contractor, to the extent 
that such employee or contractor of NASA or employee of such 
contractor prepares, d i s sem i n a te s , or provides access to any 
information pursuant to his employment or contract with NASA, or his 
employment with such contractor. 

Availabil ity Notice . Request for copies of this report should be 
referred to: 

National Aeronautics and Space Administration 
Scientific and Technical Information Facility 
P.0. Box 33 

College Park, Maryland 20740 

Final Disposition . After this document has served its purpose, 
it may be destroyed. Please do not return it to the Environmental 
Research Institute of Michigan. 


TECHNICAL REPORT STANDARD TITLE PAGE 


1. Report No. 2. Government Accession No. 

1 64000-1 9-F 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Study of Spectral/Radiometric Characteristics 
of the Thematic Mapper for Land Use Applications 

5. Report Date 

October 1985 

6. Performing Organization Code 

7. Author(s) 

William A. Mali la and Michael D. Metzler 

8. Performing Organization Report No. 

164000-1 9-F 

9. Performing Organization Name and Address 

Environmental Research Institute of Michigan 

P.0. Box 8618 

Ann Arbor, MI 48107 

10. Work Unit No. 

11. Contract or Grant No. 

NAS5-27346 

13. Type of Report and Period Covered 

Final Report 
21 Sep 1982 - 
20 Oct 1985 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, MD 20771 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Mr. Harold Oseroff (Code 902) served as Technical Officer and 
Mssrs. Brian Markham and James Irons (Code 923) served as Science 
Representatives for NASA/GSFC. 


16. Abstract 

This is a summary report of an investigation conducted in support of 
the Landsat 4/5 Image Data Quality Analysis (LIDQA) Program. Results 
of engineering analyses of radiometric, spatial, spectral, and geometric 
properties of the Thematic Mapper systems are summarized; major emphasis 
was placed on the radiometric analysis. Details of the analyses are 
presented in appendices which contain three of the eight technical papers 
which were produced during this investigation; these three, together, 
describe the major activities and results of the investigation. 


17. Key Words 

Radiometric Calibration Imaging 

Landsat-4, Landsat-5 Noise 

Thematic Mapper 
Data Quality 

Mul ti spectral Scanner Subsystem 
Information Content 

18. Distribution Statement 

Initial distribution is listed at 
the end of this document. 

19. Security Classif. (of this report) 

UNCLASSIFIED 

20. Security Classif. (of this page) 

UNCLASSIFIED 

21. No. of Pages 

49 + iii 

22. Price 


INFRARED AND OPTICS DIVISION 



Report No. 164000-19-F 


Final Report 

21 September 1982 — 20 June 1985 


on 

Study of Spectral/Radiometric Characteristics 
of the Thematic Mapper for Land Use Applications 

under 

Contract NAS5-27346 


with 


NASA Goddard Space Flight Center 
Greenbelt Road 
Greenbelt, Maryland 20771 


Submitted by 


Environmental Research Institute of Michigan 
P.O. Box 8618 
Ann Arbor, Michigan 48107 


Prepared 


William A. Malila 
Principal Investigator 



Co-Investigator 

Approved by: 

Robert Horvath 

Manager, Information Processing 
Department 



September 1985 

Original' photography may be purchase^ 
-from EROS Data Center 
Sioux Falls,. SD 5719a>^ 


INFRARED AND OPTICS DIVISION 


^ERJM 


Table of Contents 

1. OBJECTIVES 1 

2. ACCOMPLISHMENTS, CONCLUSIONS, AND RECOMMENDATIONS 1 

3. PUBLICATIONS AND PRESENTATIONS 4 

APPENDIX A: Characterization of Landsat-4 MSS and TM Digital Data. 5 

APPENDIX B: Characterization and Comparison of Landsat-4 and 

Landsat-5 Thematic Mapper Data 21 

APPENDIX C: Comparison of the Information Contents of Landsat 

TM and MSS Data 39 


i i i 


'/ 



N 86-16691 

INFRARED AND OPTICS DIVISION 


Final Report 

STUDY OF SPECTRAL/RADIOMETRIC CHARACTERISTICS 
OF THE THEMATIC MAPPER FOR LAND USE APPLICATIONS 


1. OBJECTIVES 


The objectives of this investigation were to quantify the performance of the 
Thematic Mapper, as manifested by the quality of its image data, in order to 
suggest improvements in data production and to assess the effects of the data 
quality on its utility for land resources applications. The effort was to include 
analyses of radiometric, spatial, spectral, and geometric effects, with primary 
emphasis on radiometric effects. 

This effort was part of the Landsat 4/5 Image Data Quality Analysis (LIDQA) 
program sponsored by the NASA Goddard Space Flight Center. 


2. ACCOMPLISHMENTS, CONCLUSIONS, AND RECOMMENDATIONS 


Engineering analyses were conducted to assess image data quality of the 
Landsat-4 and -5 Thematic Mappers (TMs). Macroscopic studies addressed trends 
and characteristics of full frames of data, while “microscopic” studies assessed 
differences between individual detector responses. Raw data, radiometricallv 
corrected (Type A) data, and fully corrected (Type P) data were analyzed, as well 
as special calibration data. In addition, coincident data from the TM and the 
Landsat Multispectral Scanner Subsystem (MSS) were compared. A summary of 
these analvses is given below and many of the details are presented in Appendices 
A-C. 

1. Landsat-4 and -5 Thematic Mapper (TM) data were generally found to be of 
high quality and exhibited several improvements over corresponding 
Multispectral Scanner Subsystem (MSS) data. 

a. The spatial resolution of the TM is substantially better than that of the 
MSS (nominally 30 m vs 82 m) and provided images with much greater 
spatial detail. 

b. The radiometric quantization of TM is finer (256 levels vs 64 levels at 
the spacecraft); however, part of this gain was used for increased 
dynamic range in channels rather than all being used for improved 
radiance resolution. 

c. The greater number of spectral bands of TM (seven vs four) allows 
sampling of additional spectral regions and provides increased information 
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about the scene. 

2. Detailed analyses were conducted which discovered several noticeable but 

relatively low-amplitude radiometric data anomalies or artifacts and addressed 
methods for removing them or reducing their effects: 

a. Signal Droop/Rise Along Scan Lines 

i) The mean signal level of each of the Primary Focal Plane bands 
(Bands 1-4) was seen to decay or “droop” across the scan during 
daylight data acquisition, and to “rise” with time during nighttime 
data collection. 

ii) The magnitude of the effect was found to be relatively small 
(maximum effect of less than 4% of mean signal level), but it does 
produce image banding between scan directions near frame edges. 

iii) An exponential decay model with a time constant equivalent to 
approximately 1000 pixels was found to fit the observed 
“droop/rise”; it could be used to correct data for this effect, but the 
correction would have to be computed separately for each pixel. 

b. Scan-correlated Level Shifts 

i) The mean output signal level of most detectors was found to shift 
up or down with a fixed amplitude but at essentially random 
intervals. These “scan-correlated level shifts” fell into two patterns 
for Landsat 4, named Type 4-1 and Type 4-7 after the spectral 
bands most strongly displaying each effect, and one pattern for 
Landsat 5, Type 5-3. All detectors belonging to a particular 
pattern shifted in unison, although some might shift up while 
others shifted down. 

ii) Essentially all detectors were affected, but with different 
amplitudes. The phase and amplitude relationships were found to 
be very stable, with magnitudes being as great as two video 
quantum levels (DN) in the worst case, and less than 0.5 DN in 
most other cases. 

iii) Procedures for removing the level-shift artifact from the image data 
were developed and demonstrated; they use 20-48 samples of 
calibration data which are not generally available to Landsat data 
users. Therefore corrections for this artifact should be incorporated 
in the TIPS (Thematic Mapper Image Processing System) 
processor. 

c. Quantization Step Size 

i) The quantization step sizes in the .AD converters on the Thematic 
Mappers (both Landsat 4 and Landsat 5) were unequal in size, 
resulting in uneven distribution of data within the signal space 
spanned. 

ii) Errors introduced by this artifact are small, typically one quantum 
level or less. 

iii) The resampling process used to produce fully corrected (Type P) 
data both masks the presence of this artifact and precludes a 
correction for it. 

d. Histogram Gaps (Empty Bins) 

i) The application of gains greater than unity to the raw signals 

during radiometric correction resulted in occasional “empty bins'* in 
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the individual detector histograms. 

ii) Although unappealing from an aesthetic sense, the result is a 
natural consequence of stretching the quantized data. 

iii) No correction is possible; evidence of the effect is masked in fully 
corrected (Type P) data by the detector-equalization and resampling 
processes. 

e. Coherent Noise 

i) Very low amplitude (<0.25 DN) coherent noise was observed at a 
frequency of approximately 400 Hz. 

ii) Although the source of this noise was not definitely determined, its 
effect was judged to be minimal. 

f. Thermal (Band 6) data were found to have a full-range temperature 
span of 200°K to 340°K instead of the advertised 260°K to 320°K; the 
correct range is indicated in the calibration parameters provided in the 
digital tape headers. 

3. Radiometric comparisons were made using data from simultaneous acquisitions. 

a. Simultaneous coverage by Landsat 4 and Landsat 5 Thematic Mappers 
occurred during orbit adjust maneuvers shortly after the launch of 
Landsat 5 and allowed direct radiometric comparison of the two sensors 
to be made. 

i) Relationships were established between data from the two sensors; 
gain differences as large as 14 percent were noted in calibrated 
data. 

ii) Low-level clipping was observed in the radiometrically corrected 
Landsat-5 data for Bands 5 and 7, but not in raw data. 

iii) A continued monitoring of system performance and data calibration 
would be valuable to many users of Landsat data. 

b. Simultaneous Landsat-4 TM and MSS coverages allow comparisons of 
the radiometry from those two sensors as well; one such scene pair was 
analyzed. 

i) High degrees of correlation were found; relationships to convert 
single spectral bands of TM data into MSS-equivalent data were 
developed. 

ii) Spectrally transformed features (Tasselled-Cap Brightness and 
Greenness) were compared; Greenness values were almost perfectly 
correlated, but Brightness values exhibited some differences 
(correlation coefficient of 0.75). 

iii) TM data were observed to have higher spectral dimensionality than 
the MSS data. 

c. Information theoretic measurement techniques were developed and applied 
to simultaneous TM and MSS data. Analyses of the data volumes 
realized by each sensor relative to the information potential of the 
sensors were made, and the greater information content of TM data 
from an agricultural scene was quantified. 
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3. PUBLICATIONS AND PRESENTATIONS 


In addition to the ten Quarterly Status and Technical Progress Reports, and 
presentations at Investigator Workshops and Symposia, eight technical papers 
dealing with efforts under this contract have been published in journals and 
symposium proceedings. They are listed below. Copies of Papers 4, 7, and 8 are 
included as Appendices A, B, and C, respectively. These three papers, together, 
describe the major activities and results of the investigation. 

1. Metzler, M.D., and W.A. Malila, “Scan-angle and Detector Effects in Thematic 
Mapper Radiometry,” in Proc. Landsat— 4 Early Results Symp., NASA Goddard 
Space Flight Center, Greenbelt, MD, February 1983. 

2. Metzler, M.D., and W.A. Malila, “Radiometric Characterization of Thematic 
Mapper Full-Frame Imagery/' in Proc. SPSE/APS Conf. on Techniques for 
Extraction of Information from Remotely Sensed Images , Rochester, NY, August 
1983. 
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Falls, SD, Octoberr 1983. 
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Data from Landsat MSS and TM,” in Proc. 18th International Symp. on 
Remote Sensing of Environment , Paris, France, October 1984. 
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pp. 1449-1457. September 1985. 
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Abstract- Engineering analyses conducted to assess image data quality 
are described and results are presented for the Landsat-4 Multispec- 
tral Scanner (MSS) and Thematic Mapper (TM). Also, coincident data 
from the two sensors are compared. Macroscopic studies addressed 
trends and characteristics of full frames of data, while “microscopic" 
studies assessed differences between individual detector responses. Raw 
data, radiometric ally corrected data, and fully corrected data were 
analyzed, as well as special calibration data. 

The Landsat-4 MSS was found to produce data of high quality, 
comparable to previous Landsats, except for a low-level coherent 
noise effect which is unique to the current sensor. Radiometric rela- 
tionships between Landsat-2 and -3 MSS’s and the Landsat-4 MSS 
were established through empirical analysis of simultaneously acquired 
data. 

The TM was found to produce image data of very good spatial reso- 
lution and overall good radiometric data quality, showing improvements 
over MSS. Radiometric equalization of detector responses was found 
to be close to the theoretical limit of quantization error, except for two 
relatively low-amplitude artifacts. One is a difference in response that 
depends on the direction of mirror scan. This produces scan-angle 
effects superimposed on scene-related effects. The second is a tendency 
for level shifts to occur between mirror scans at random times but with 
correlations between detector responses. Two forms or patterns of 
level shift were identified, corresponding to four system noise states. 
Preliminary correction models and/or procedures have been developed 
and recommended for further evaluation. 

In a comparison of coincident coverages of one agricultural scene, 
high degrees of correlation were found between MSS signals and those 
from corresponding TM bands. Very high correlation also was observed 
between spectrally transformed Greenness variables for the two sensors. 
Lower correlation was observed for Brightness. 

Keywords - Remote Sensing, Landsat-4, Multispectral Scanner, 
Thematic Mapper, Digital Image Processing. Radiometric Calibration. 
Coherent Noise, Scan Effects, Level Shifts, Detector Equalization, Spa- 
tial Resolution. 


I. Introduction 

T HE JULY 1982 launch of Landsat-4 and its two sensors 
represented both a continuation and a new beginning 
in the remote sensing of Earth resources [ I ] . A second decade 
of satellite imaging with the Multispectral Scanner (MSS) is 
underway. Its 80-m spatial resolution is fine enough to resolve 
many natural features and land-use details in both rural and 
urban settings. Coupled with the spectral information of its 
four bands, many applications have evolved and developed 
during the past decade. A new era of sensing with refined 

Manuscript received December 2, 1983; revised December 27, 1983. 
This work was supported by the National Aeronautics and Space Ad- 
ministration (NASA), Goddard Space Flight Center, Greenbclt, MD 
under Contracts NAS5-27254 and NAS5-27346. 

The authors are with the Environmental Research Institute of Michi- 
gan, Ann Arbor, MI 48107. 


spatial resolution (30 m) and expanded spectral coverage (7 
bands) has begun with the orbiting of the Thematic Mapper 
(TM) [2] , [3] . Although telemetry problems have restricted 
the acquisition of data by the Landsat-4 TM, sufficient data 
have been analyzed to excite investigators about its increased 
information potential and evoke anticipation for the forth- 
coming launch of Landsat-D’ and the second TM sensor. 

This paper describes results from engineering studies of the 
characteristics of digital image data from the two Landsat-4 
sensors [4] -[7]. These studies are part of the Landsat-4 
Image Data Quality Analysis program (LIDQA) which was 
organized and sponsored by the National Aeronautics and 
Space Administration (NASA), Goddard Space Flight Center 
(GSFC), Greenbelt, MD [8] . Control of the MSS sensor and 
system was transferred from NASA to the National Oceanic 
and Atmospherics Administration (NOAA) in January 1983, 
while the TM system is planned to remain on experimental 
status until a later date. The studies described here were fo- 
cussed on the radiometric characteristics of the MSS and TM, 
with limited consideration of spectral and spatial properties. 
Studies such as these are desirable so that subsequent investi- 
gators and users of data will have verification of the system’s 
performance and knowledge of any artifacts, without the 
necessity of each conducting his or her own detailed analysis. 
Also, NASA uses the results of internal and external LIDQA 
investigators to make improvements in product generation 
procedures. 

We believe that continuity of an operational remote sensing 
system, like the Landsat MSS, has at least two aspects: a) 
continued availability of high-quality data and b) consistent 
interpretability of signal amplitudes and spectral features ex- 
tracted from those data. 

Both sensors scan linear arrays of detectors perpendicular 
to their ground track in order to achieve coverage of the scene 
beneath. If detectors within a band are not perfectly calibrated 
with respect to each other (i.e., equalized), image banding or 
striping results and extra variability or noise is introduced into 
the statistical descriptions (e.g., signatures) of signals from 
specific scene classes. Similarly, if the noise associated with 
each detector channel were to significantly exceed design 
specifications, the utility of the data collected by the system 
would be diminished. Also, reduced dynamic range could 
degrade data utility. 

In addition, the sensor response should be consistent and 
preferably constant as a function of scan angle, position along 
track, and days in orbit, to permit unambiguous and consistent 
information extraction and interpretation. For MSS, this ex- 
tends to a need either to match the radiometric calibration of 
the preceding Landsat MSS sensors or to establish relation- 
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ships for appropriately adjusting coefficients and thresholds in 
affected algorithms and procedures that were developed for 
earlier data. For TM, which employs bidirectional scanning, 
one would not want response to depend on scan direction. 

Continuity in spectral response is important to MSS. It also 
is desirable to compare the spectral information content of 
TM to that of MSS in order to assess both the extent to which 
the same information can be extracted from TM and any in- 
crease in information that may be available because of the 
additional spectral bands in TM. 

II. Approach and Method 

Both macroscopic and microscopic analyses were conducted 
of digital computer-compatible-tape (CCT) data from both 
sensors. Data from several stages of the ground data processing 
procedure were analyzed, including raw data, radiometrically 
corrected data, fully corrected data, and calibration data, as 
available. 

The macroscopic analyses were designed to assess trends 
within spectral bands over the full 185 X 185-km frames of 
data. First, histograms and statistical measures (e.g., means 
and standard deviations) of signal amplitudes were computed 
for each band. Where available, radiometrically corrected and 
fully corrected duta were examined and compared. Selected 
scatterplots were made of data in pairs of bands and signal 
correlations computed. MSS responses from selected scene 
classes also were compared qualitatively with expectations 
based on our experience with previous Landsats. 

Second, scan-angle effects were investigated by averaging 
columns of pixels down the entire frame, plotting the results 
(amplitude versus pixel number or scan angle) for each band, 
and analyzing these average scan lines. Third, to explore the 
other image dimension, down-track profiles (amplitude versus 
scan line number) were computed, one for each band, by aver- 
aging all pixels (or a systematic sample) on each line to obtain 
a single average value for each scan line. Fourth, to study any 
scan direction effects in TM, average scan lines were computed 
for each band, separately for each of the two scan directions. 
Computer plots and listings of the computed average scan lines 
and down-track profiles were produced and analyzed. 

Finally, to compare radiometric calibration between the 
Landsat-4 MSS and Landsats-2 and -3, opportunities for 
simultaneous acquisitions were identified and a few pairs 
of scenes were acquired for analysis. Simultaneous acquisition 
was possible on selected dates at selected locations due to the 
16-day cycle of Landsat-4 and the 18-day cycle of Landsats-2 
and -3. 

The microscopic analyses were designed to assess differences 
between individual detector responses before and/or after 
equalization within each spectral band (6 detectors for each 
MSS band, 16 for each TM reflective band, and 4 for the TM 
thermal band). They also were designed to examine noise 
characteristics by detector, both down-track and along the 
scan line. Many of the analyses conducted on spectral band 
values were repeated for each individual detector to provide 
data for comparison. Detector histograms and related statistics, 
together with the down-track profiles, permitted assessment of 
the ground systems’ success in equalizing detector responses. 


A unique aspect was our analysis of nighttime TM data from 
the reflective channels. Fourier analysis of down track profiles 
also was used to assess the success of equalization and to de- 
tect the presence of periodic banding and striping effects. A 
Fast Fourier Transform (FFT) technique adapted from [9] 
was used. 

A coherent noise effect was visually apparent in MSS data. 
Fourier analysis of individual scan lines, selected throughout 
the frame, was accomplished using the FFT procedure. The 
frequencies of anomalously large components were identified 
and their amplitudes measured from the FFT analysis output. 

Geometric and spatial characteristics also were considered. 
The pixel lengths of MSS scan lines were examined and com- 
pared in radiometrically corrected data (CCT-AT) from an 
Eastern U.S. scene acquired when TM was operating and a 
Western U.S. scene acquired when TM was off. For comparison 
of MSS and TM and for illustration of spatial-resolution ef- 
fects on one’s ability to resolve fine detail, an appropriate 
scene (North Carolina, Path 14/Row 36) for which simul- 
taneous MSS and TM coverage existed was selected. Scan line 
traces (amplitude versus pixel number) were plotted in a 3-D 
projection to permit comparison and to illustrate differences. 
Quantitative comparison was based on computed signal vari- 
ances within individual fields. 

Finally, a limited comparison of spectral responses from 
MSS and TM was made for the same simultaneously acquired 
scene. Spatial registration was accomplished by means of pixel 
shifts based on control points identified in the two images; an 
accuracy of approximately one pixel was achieved. Groups of 
four TM pixels were averaged for comparison with their asso- 
ciated MSS pixels. Correlations were computed between MSS 
and comparable TM bands. Also, correlations were computed 
between spectrally transformed features, namely Brightness 
and Greenness, defined for the two sensors [10]. The MSS 
transformation was equivalent to the Tasseled-Cap Transfor- 
mation of Kauth and Thomas [11], while the TM transforma- 
tion is an analogous one developed by Crist and Cicone [12] . 
The reader is referred to [12] for a more detailed discussion 
of the spectral characteristics and spectral information content 
ofTM. 


III. Data Sets 

Table I identifies the 9 MSS frames which were analyzed 
part of this study ; they were acquired during the first four 
months of Landsat-4 operation. The 6 TM frames analyzed 
are listed in Table II; they represent the first six months of op- 
eration. Radiometrically corrected (CCT-AT) digital data tapes 
were the primary data sources for our MSS studies, supple- 
mented by fully corrected imagery. Raw data (CCT-BT) and 
radiometrically corrected data (CCT-AT) were the primary in- 
puts for our TM studies, with added analyses of fully (geomet- 
rically as well as radiometrically) corrected data (CCT-PT) and 
calibration data (CCT-ADDS). The CCT-BT and CCT-ADDS 
are special engineering products not generally available to users. 
The MSS data were products of the Multispectral Scanner 
Image Processing System (MIPS), while the TM data were pro- 
cessed through the interim processing system (SCROUNGE) 
established at NASA/GSFC. 
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TABLE I 

Mss Analysis Data Set 


Geographic 

Location 

Landsat 

Path/Row 

Date 


Frame 

S. Carolina 
(Inland) 

4 

17/36 

29 Sep 

82 

40075-15271 

California 
( Imperial 
Valley) 

4 

39/37 

23 Sep 

82 

40069-17433 

N. Carolina 
(coast) 

4 

14/36 

24 Sep 82 

40070-15084 

N. Carolina 
(coast) 

3 

15/35 

24 Sep 

82 

31664-15070 

- 

4 

14/35 

24 Sep 

82 

40070-15081 

New England 

3 

14/30 

22 Dec 

82 

31753-14591 

• 

4 

13/30 

22 Dec 

82 

40159-15010 

New Mexico 

2 

34/37 

9 Nov 

82 

22848-16571 

- 

4 

32/37 

9 Nov 

82 

40116-17005 


TABLE II 

TM Analysis Data Set 


Geographic 

Location 

Path/Row 

Oate 

Frame Number 

Scene Type 

Arkansas 

23/35 

22 Aug 82 

40037-16034 

Agriculture, some 
clouds, water 

Iowa 

27/30 

3 

Sep 82 

40049-16262 

Agriculture, clear 
skies 

North 

Carolina 

14/36 

24 

Sep 82 

40070-15084 

Agriculture, much 
cloud, water 

Cape Cod 

11/31 

8 

Dec 82 

40145-14492 

Low radiance, clouds, 
much water 

Georgia 

116/207 

24 

Dec 82 

40161-02481 

Nighttime 

Grand 

Bahamas 

14/42 

14 

Jan 83 

40182-15125 

water, clouds 


IV. MSS Analysis Results [4], [5] 

A. Landsat-4 MSS Data Quality 

1) Qualitative Comparisons with Data from Preceding 
Landsats: Scatterplots and statistics derived from data from 
the first two scenes listed in Table I were examined, in both 
comprehensive samples from the entire frames and samples 
from distinct scene classes. Fig. 1 displays the MSS3 -versus* 
MSS2 data dispersion of the comprehensive sample from the 
Imperial Valley scene; spectral regions occupied by four scene 
classes are outlined on the figure. Water is represented by low 
signals in all bands, while sand is bright in all bands. Clouds 
were so bright as to cause signal saturation (Level 127) at the 
prevailing solar illumination angle of 47° (elevation). Also, 
agricultural data can be seen to fill out a roughly triangular 
shape (a rotated “tasseled cap” shape). These observations 
are generally true also for data we had examined previously 
from Landsats-1-3. 

In addition, we compared data from the simultaneously 
acquired Landsat-3 and Landsat-4 data pair over the coast 
of North Carolina. Image displays and general statistics showed 
good correspondence. 

2) Detector Equalization and Quantization Effects: Visual 
examination of photographic images and other displays of 
MSS data revealed evidence of stripes and bands caused by 
incomplete detector equalization, in areas of nearly uniform 
signals (such as water bodies and barren plains). Signal means 


computed for each detector in five such areas of the Imperial 
Valley scene are presented in Table III. Among the six detector 
means for each band, the standard deviation was generally less 
than 0.3 counts, while 0.88 counts was the largest difference. 
A more comprehensive estimate of residual banding and striping 
throughout the first two MSS frames listed in Table I came 
from the FFT analysis of down-track profiles. Evidence of 
radiometric disturbances, which did not appear to be asso- 
ciated with ground phenomena, were observed at spatial wave- 
lengths of 2, 3, and 6 pixels. The largest calculated average de- 
viation from the mean of the 6 detectors was 0.27 ± 0.05 count 
for the South Carolina frame. This result is in good agreement 
with the Imperial Valley values in Table III. 

The detector equalization algorithm used by the production 
system for full-frame MSS data sets [13] was examined to see 
if any further improvement seemed possible. The algorithm 
carries out signal decompression (Bands 1-3 only), calibration 
wedge normalization, and detector histogram equalization 
using a single transformation of received satellite samples to 
determine final count levels. A separate transformation is 
computed for 600-scan-line segments within the 2400-scan- 
line frame, and then is interpolated to each 200-scan -line 
sub-segment. The result is a pattern of full count-level bins 
with occasional empty ones, a pattern which differs among 
detectors due to both their differing responses and the correc- 
tion process. Examination of this pattern for a series of 200- 
scan-line sub-segments showed the expected gradual change 
from one segment to the next, and indicated that the algorithm 
worked as intended. 

A very important limit on accuracy of this equalization algo- 
rithm is not one of measuring differences or determining cor- 
rections, but rather one of translating corrected values to non- 
fractional count levels. The theoretical rms error of this 
truncation process is 0.29 [14] , which is similar to the error 
magnitudes reported above. That is, the current correction 
algorithm cannot be improved much if any without dealing 
with the truncation error problem. This error could be re- 
duced somewhat (per radiance unit) by using more numeric 
precision, for example by spreading the signal values over the 
range 0-255 instead of the current 0-127. However, the 
analog-to-digital conversion process in the satellite and its 
similar assignment to nonfractional values during its six-bit 
quantization can cause minor but uncorrectable quantization 
effects which will be most noticeable in relatively uniform 
areas. 

3) Coherent Noise: A careful visual examination of imagery 
from Landsat-4 MSS revealed a diagonal striping pattern, 
particularly in areas of nearly uniform radiance. Our FFT 
analysis of this noise pattern on selected scan lines distributed 
throughout the first two Landsat frames of Table I revealed 
that the noise existed at a dominant spatial, wavelength of 
about 3.6 pixels along scan, as illustrated by Fig. 2, and was 
present in all bands. The frequency corresponding to the 
observed dominant wavelength is approximately 28 kHz. The 
wavelength was slightly different (3.57 compared to 3.59) for 
the two frames, but was consistent (± .001) within each frame. 
In one frame (South Carolina), two other wavelengths appeared 
somewhat consistently, namely the wavelengths of 2.02 and 
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Fig. 1. MSS spectral distribution of scene classes from the Imperial Valley, CA scene. 


TABLE III 

Detector Averages for 5 Areas in One LANDSAT-4 MSS Frame 


(Imperial Valley. CA) 


Detector 


1 

2 

3 

4 


1 


23.01 

21.08 

36.80 

14.44 


2 


23.29 

21.00 

36.91 

14.41 


3 


22.94 

21.10 

36.94 

14.53 


4 


23.23 

21.04 

37.31 

14.62 

Highlands 

5 


23.25 

20.97 

37.18 

14.43 


6 


23.31 

21.03 

37.15 

14.70 


Std. Dev. 

* 

.16 

.05 

.19 

.12 


1 


33.08 

35.51 

53.86 

19.03 


2 


33.38 

35.78 

54.18 

19.09 


3 


33.41 

36.06 

54.32 

18.99 


4 


33.32 

36.10 

54.22 

18.80 

Aqrlcul tural 

5 


33.31 

35.86 

53.60 

18.81 


6 


33.08 

35.69 

53.75 

18.91 


Std. Dev. 

• 

.15 

.22 

.29 

.12 


1 


62.23 

76.10 

84.14 

25.73 


2 


62.11 

76.76 

84.25 

25.62 


3 


62.10 

76.08 

84.54 

25.84 


4 


61.42 

76.01 

84.04 

25.89 

Sand 

5 


61.78 

75.88 

83.81 

25.84 


6 


61.99 

76.23 

83.88 

25.63 


Std. Dev. 

• 

.30 

.31 

.27 

.12 


1 


20.79 

13.65 

10.09 

2.43 


2 


21.02 

13.69 

9.88 

2.35 


3 


20.64 

13.70 

9.91 

2.21 


4 


20.73 

13.52 

9.59 

2.27 

water 

5 


20.26 

13.37 

9.84 

2.07 


6 


20.81 

13.90 

10.22 

2.30 


Std. Dev. 

* 

.25 

.18 

.22 

.12 


1 


51.62 

61.55 

69.43 

22.10 


2 


51.77 

61.82 

69.17 

22.27 


3 


51.58 

61.58 

69.30 

22.44 

Bright Shallow 

4 


52.16 

61.68 

69.52 

22.53 

Valley 

5 


51.28 

62.07 

69.31 

21.97 


6 


51.67 

61.79 

69.56 

22.05 


Std. Dev. 

» 

.29 

.19 

.15 

.22 



4.63 pixels, but these and other smaller peaks in both frames 
were not analyzed in detail. The magnitude of the sine-wave 
component at the primary peak was computed and tabulated 
for the group of scan lines that was processed. This informa- 
tion, presented in Table IV, shows that the rms magnitude of 
the noise did not exceed one count in the worst band (MSS1). 

NASA scientists now suspect that this effect may be due to 
a coherent noise source in a power supply. Based on the mag- 


nitudes described above, we consider this effect to be a rela- 
tively minor source of error, in spite of its nuisance effect on 
image appearance and its possible effects on classification re- 
sults for spectrally similar classes. Nevertheless, it would be 
desirable to eliminate this noise source from future sensors. 

4) Geometric Distortion: In carrying out the MSS radio- 
metric analysis described above, a significant geometric dis- 
tortion was noted in some CCT-AT tapes used in the analysis. 
This distortion was characterized by adjacent groups of 6 scan 
lines (each group comprising one mirror sweep) being dis- 
placed relative to one another, as illustrated in Fig. 3. The 
magnitude of the displacement varied up to 6 pixels and de- 
creased to zero at the left side of the image, indicating that the 
distortion is a line length variation. We found that this distor- 
tion was absent or much reduced in one scene (Imperial Valley) 
taken when the Thematic Mapper sensor was not operating. 

For further analysis of the scenes acquired with the TM 
operating, we implemented a linear line stretch correction, 
based on the position of end-of-scan codes in each scan line, 
which we understand is equivalent to the correction used by 
the EROS Data Center (EDC) to produce the geometrically 
corrected Landsat-4 CCT-PT tapes. This correction substan- 
tially improved the geometric correspondence between paired 
Landsat-3 and -4 datasets. As Fig. 3 shows, this correction re- 
stored the river .(vertical band of dark symbols) to its proper 
contiguous course. Thus fully corrected data from EDC should 
exhibit little if any residual effect, but users of CCT-AT tapes 
should expect to find displacements. 

B . Between-Satellite Radiometric Calibration of MSS 

1) Matched Datasets: Three loci of frame centers in the 
U.S. were identified along which ground paths of Landsats-3 
and -4 coincided sufficiently to permit similar view directions 
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WAVELENGTH (Pixels) 

Fig. 2. Along-scan Fourier transform showing MSS coherent noise peak at wavelength 3.6 (South Carolina Scan Line 

2117. MSS3). 


TABLE IV 


Average Magnitude of MSS Coherent Noise Effect 


Band 

Maqmtude (counts) 

1 

0.75 ♦ .11 

2 

0.52 ♦ .09 

3 

0.56 ♦ .10 

4 

0.50 ♦ .10* 


•This Band 4 value has been scaled to reflect the calibration of 
data oroduced on or after 20 October 1982. 


to scene locations. Fig. 4 shows one along the East Coast, 
one in the central U.S. and the other near the West Coast. The 
two outside loci would allow repeated coincident coverages 
every 144 days (144 being the least common multiple of 18 
and 16), but coverages no closer than one day apart were 
possible for the central locus. Several of the identified poten- 
tial coincident-coverage datasets were requested and approved 
lor collection. After losses due to cloud cover and acquisition 
problems were incurred, collection of coincident Landsat-3 
and -4 data was successful on 5 dates in the U.S. Further- 
more, although Landsat-2 has not been collecting MSS data 
operationally since February 1982, it was brought to our at- 
tention that a coincident Landsat-2/Landsat-4 dataset was 
collected in November 1982 and we were provided a copy 
tor analysis [15]. The coincident-coverage datasets we ana- 
lyzed are listed in Table I ; regression relationships are presented 
in the next section for the New England pair for Landsats-3 
and -4 and for the New Mexico pair for Landsats-2 and -4. 

21 Empirical Regression Relationships: From the New 
England Landsat-3/4 data set, 60 relatively large and uni- 
form polygonal areas were carefully identified and located in 
both frames. These areas were selected to provide as wide a 
spectral range as possible. Table V presents empirically derived 



After Correction 

Fig. 3. Line-length variation in Landsat-4 MSS CCT-AT tapes (South 

Carolina). 


linear regression coefficients for the following equation: 

[Landsat-3 Signal) = Gain * [Landsat-4 Signal] + Offset. 

Note that very high R 2 values and small standard errors on 
the order of one-half count were achieved. Also included in 
the table are coefficients based on published pre-launch cali- 
bration values. These pre-launch gains are seen to be approxi- 
mately 10 percent lower than the empirical values, for data 
processed by the EROS Data Center (EDC) between October 
20, 1982 and April 1, 1983; small adjustments in EDC coef- 
ficients were made on April 1, 1983 [16]. 
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Fig. 4. Loci of frame centers with potential coincident coverage by Landsats-3 and -4. 


TABLE V 

Empirical LANDSAT-4 to -3 Calibration Coefficients 


Scene 

Band 

Gain 

Offset 

SE* 

d 

New England 

1 

1.018 

-1.614 

0.52 

0.996 


2 

1.112 

0.018 

0.45 

0.999 


3 

0.9096 

-0.463 

0.53 

0.999 


4 

1.148 

-0.421 

0.53 

0.998 

Prelaunch 

1 

0.894 

-1.00 



Calibration** 







2 

1.000 

0.722 




3 

0.863 

0.870 




4 

1.026 

2.34 




Notes: 

*SE * standard error In slqnal counts 

••This Information represents carefully conducted ground measure- 
ments taken before launch of each satellite, along with 
aoorooriate scaling of the data to account for ground 
processing algorithm coef flcients. 

•••The values In this table assume EROS Oata Center ( EOC ) processing 
of the data between 20 October 1983 and 1 April 1983; EOC sources 
should be consulted for coefficient adjustments on and after 1 
Aorll 1983. 

Relationship: 

Landsat 3 • Gain x (Landsat 4) ♦ Offset 

To contrast the effects of using various forms of data cor- 
rection, potential errors were calculated for dark, interme- 
diate, and bright areas. Between 2 and 12 counts of error 
would occur, using an Euclidean count difference measure, if 
Landsat-4 data were to be used directly as a substitute for 
Landsat-3 data. If pre-launch measurements of calibration 
were used to correct the difference, the error would still fall 
between 3 and 9 counts, an insubstantial improvement. Thus 
use of the empirical relationships in Table V is recommended, 
with the caution that any additional changes in ground pro- 
cessing coefficients should be taken into account. 

Table VI presents an analogous set of coefficients for Land- 
sats-2 and -4, based on 15 polygonal areas. Note that the stan- 
dard errors are approximately twice as large as for Landsats-3 
and -4, indicating less accuracy. Besides being fewer in num- 
ber, the polygons used for Landsat-2/4 were more difficult to 
locate, and there were patterns of haze variation that could 


TABLE VI 

Empirical LANDSAT-4 to -2 Calibration Coefficients 


Scene 

Band 

Gain 

Offset 

SE 

R- 

New Mexico 

1 

0.8766 

-0.592 

0.96 

0.973 


2 

1.0124 

2.187 

1.12 

0.989 


3 

0.8796 

1.699 

0.97 

0.995 


4 

1.1002 

-2.629 

0.59 

0.997 

Prelaunch 

1 

0.894 

-2.988 



Cal ibration 

2 

1.035 

-1.494 




3 

0.863 

-1.740 




4 

1.026 

-0.334 




Relationship: 

(Landsat 2) * Gain x (Landsat 4) + Offset 
See Notes under Table V. 

possibly have caused some differences when coupled with a 
2.9° difference in view angle and a 4-min difference in acquisi- 
tion time. 

V. Comparison of MSS and TM 
A. General Characteristics and Specifications 
Table VII summarizes the major differences between MSS 
and TM specifications, indicating ways in which TM (and/or 
Landsat-4) represents extensions of the spatial, spectral, tem- 
poral, and radiometric characteristics of the MSS of Landsats 
1-3. Specific details are readily available from [l]-[3] or 
related papers in this issue. 

These improved specifications offer the potential for greater 
information gathering capacity, given appropriate information 
extraction techniques. LIDQA engineering studies are directed 
at assessing the extent to which particular performance goals 
have been met and searching out any unwanted system charac- 
teristics that might interfere with extraction of the information. 
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TABLE VII 

Extensions Offered by TM/ LANDSAT-4 Beyond the 
Characteristics of MSS/LANDSATs 1-3 

1) SPATIAL SPECIFICATIONS 

a) 30m vs. 80 m resolution 

t>) I mo roved platform stability (Landsat 4 vs. Landsats 1-3) 

2) SPECTRAL SPECIFICATIONS 

a) Narrower soectral bandwidths 

b) Additional soectral regions (Blue-green, intermediate IR, 

thermal IR) 

3) TEMPORAL SPECIFICATIONS (LANDSAT 4 VS. LANOSATS 1-3) 

a) Repeat cycle of 16 vs. 18 days 

b) Reoeat coverage of overlap zone of 7 or 9 days vs. next 

day 

4) RADIOMETRIC SPECIFICATIONS 

a) Imoroved slgnal-to-nolse 

b) Increased dynamic range (8 bits vs. 6 bits) 


B. Comparisons of Coincident Data 

1) Site Description: The site selected for our limited, but 
direct, comparison of MSS and TM data [10] is a private 
agricultural holding along the coast of North Carolina. This 
recently cleared area of 18 000 hectares is divided into square- 
mile blocks by improved two-lane gravel roads. Each block is 
bordered by drainage ditches and is divided into 16 strip fields 
(each 100 m wide) by a series of 3-m-wide drainage ditches. 
Fields of soybeans, corn, and fescue at various stages of de- 
velopment were present. An unusual amount (nearly 20 cm) 
of rain had fallen just prior to the overpass data Sept. 24, 1982. 
MSS2 and TM3 images of the site are presented in Fig. 5. 

2) Spatial Comparison: The images of Fig. 5 very clearly 
illustrate the improved resolving power of the IM system. Field 
edges, roads, and other features are sharp and well defined. In 
many of the large square-mile fields (those having little vegeta- 
tion present), the patterns of drainage ditches and narrow strip 
fields are clearly seen in the TM image but not in the MSS. 
Another graphical display of the differences in resolving power 
of the two systems is presented in Fig. 6. A computer plotter 
produced a horizontal line for each pixel at a height propor- 
tional to its amplitude. Lines connecting the levels of adjacent 
pixels give rise to a random staircase -like pattern for each scan 
line illustrated. Multiple lines crossing Test Region 1 (outlined 
on Fig. 5) are displayed in a 3-D perspective and clearly show 
the strip-field pattern for TM but not for MSS. Histograms 
and signal variances substantiated and quantified this difference. 

3) Spectral Comparison: Three TM spectral bands (TM2, 
TM3, and TM4) span nearly the same spectral region as the 4 
MSS bands. The spectral diversity of the test site provided a 
good data set for computing correlations between data values 
from these two sets of bands. Results are shown in Table VIII. 
Correlations of TM4 were about 0.97 and 0.99 between MSS3 
and MSS4, respectively, and much lower and of opposite sign 
for the two visible bands, as would be expected due to vegeta- 
tion characteristics. TM3 was most highly correlated with 
MSS2 (0.98), and quite well correlated with MSS1. Finally, 
TM2 was well correlated (0.91) with both MSS1 and MSS2. 
Table VIII also contains correlations between Brightness and 
Greenness variables computed for the two data sets. The 
Greennesses are very highly correlated (0.99), while the Bright- 


nesses are less well correlated (0.75), a result which warrants 
further investigation and is likely due in part to the added 
information of the intermediate IR bands in TM Brightness. 

VI. TM Analysis Results [6], [7] 

A. General Characteristics of TM Data 

1) Signal Histograms and Dynamic Range: Signal values 
over the full 0-255 signal range have been observed in all 6 
reflective bands, with clouds frequently causing saturation. 
However, in a cloud-free, primarily agricultural scene (Iowa), 
we found narrow concentrations of signals in all bands but 
TM4 and TM5 which had somewhat broader distributions 
(see Fig. 7). 

Fig. 8 presents several TM4 signal histograms to illustrate 
other TM characteristics. In parts (a), (b) and (c), responses 
from all detectors are grouped together. A ragged histogram 
shape is evident for the raw data of part (a), the result of un- 
equal step sizes in the analog-to-digital converter of the sensor. 
After radiometric correction, part (b), the levels of individual 
detectors are shifted with respect to each other and the ragged- 
ness is decreased. Finally, the cubic -convolution interpolation 
process used in geometric correction results in a smooth histo- 
gram, part (c), for the fully corrected data. Parts (d) and (e) 
are for a single detector, with part (d) again displaying the A/D 
step size variation. When the radiometric correction process 
involves a stretching of the signal range (gain > 1 .0), occasional 
gaps (empty levels) are created in the histogram as a direct 
result of the quantized, integer nature of the signal values, as 
shown in part (e). Because the detectors have slightly dif- 
ferent gains and offsets, these empty levels are not evident in 
the histogram for the entire band, part (b), although they will 
still be present for individual detectors in CCT-AT data. 

2) Detector Equalization and Down-Track Profiles: Just as 
for MSS, detector equalization within TM bands was found to 
be close to the theoretical limit of quantization. Residual dif- 
ferences in CCT-AT detector means were evident in corre- 
sponding CCT-PT scan-line means (down-track profiles) and 
in imagery. 

In addition to the residual differences in detector means, we 
discovered another down-track artifact: very low frequency, 
scan-related noise in the form of level shifts that occurred in 
the mean values of individual detector means [6] . Other in- 
vestigators also reported related observations [17] . This arti- 
fact is discussed in more detail in Section VI-C. 

3) Scan-Angle Effects: Upon examination of plots of 
average scan lines (amplitude versus scan angle (pixel number)), 
it was found that values near the Western edge of the scene 
were substantially higher than those near the Eastern edge, as 
shown in Fig. 9. The effect, more pronounced in the short 
wavelength and thermal bands, can be attributed to sun-view 
geometry and a combination of atmospheric, bidirectional- 
reflectance, and shadowing effects. Previous modeling work 
has predicted atmospheric scan angle effects similar to those 
observed here [18] , [19] . A plot of model results derived by 
us from the data set of Dave [20] is illustrated in Fig. 10. As 
the wavelength increases, atmospheric effects are reduced and 
the reflectance effects become the drivers of observed trends. 
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(b) 

Fig. 5. Comparison of MSS and TM images of agricultural scene (1 identifies the region plotted in Fig. 6). (a) MSS Image. 

(b) TM Image. 


Shading and view direction contribute to the observed thermal 
responses. 

In addition to the expected scene-related scan-angle effects 
described above, we discovered a sensor-related scan -angle ef- 
fect that is associated with the direction of mirror scan [6] . 
This artifact is discussed in more detail in Section VI-B. 

13 


B. Scan-Direction Artifact and Correction Model 
1 ) Description of Artifact: From plots of average scan lines 
for the two scan directions, we found a slight but consistent 
anomaly, the difference between West and East edge values 
was greater for forward (W-E) scans than for reverse (E-W) 
scans. Fig. 11 illustrates the results for Bands TM1 and TM6 
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Fig. 6. Geometric resolution comparison -Test Region 1: BRIGHT- 
NESS. (a) MSS-Perspective #1. (b) MSS-Perspective #2. (c) TM- 
Perspective #1. (d) TM-Perspective #2. 


TABLE VIII 

Correlation Results for MSS and TM 
(Averaged Pixels) 


a) Original 8and Values 


TM 2 

.909 

.907 

-.209 

-.356 

TM 3 

.881 

.975 

-.455 

-.574 

TM 4 

-.391 

-.586 

.972 

.986 


MSS 1 

MSS 2 

MSS 3 

MSS 4 


b) Tasseled-CaD Features 

TM Brlghtness/MSS Brightness .747 

TM Greenness/MSS Greenness .990 


for the Iowa frame, but it was observed in other bands as well. 
The previously described scan-angle effect is still dominant, 
but it appears that a systematic droop of signal values occurs 
during the active scan. This tends to increase the overall 
scan-angle effect for forward scans and reduce it for reverse 
scans for the reflective bands, being most pronounced in TM1. 
Thus 16-detector-wide swaths (approximately 17 lines on 
CCT-PT tapes and images) can differ from neighboring swaths 
by opposite amounts at the two frame edges, while little if any 
of this effect will be present at frame center. 

In the thermal band, TM6, the swaths are four detectors 
wide and, as seen in Fig. 11, the scan-direction artifact for the 
September 1982 scene was markedly different from that of 
Band TM1 and the other reflective bands. This unique thermal 
band artifact has been observed in other scenes and by other 
investigators [23] . The problem is under study by NASA. 

2) Correction Model for Band TM1: In order to charac- 

M 


CCT-AT 4-0049- 16-262 TM P027 ROJO 
BAND 1 



CCT-AT 4-0049-16-282 TM P027 ROJO 
BAND 2 



(b) 


CCT-AT 4-0049-16-262 TM P027 ROJO 
BA NO J 



Fig. 7. Full-frame histograms for TM (Frame 40049-16262). (a) Band 
1. (b) Band 2. (c) Band 3. (d) Band 4. (e) Band 5. (f) Band 6. 
(g) Band 7. 
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Fig. 7. Full-frame histograms for TM (Frame 40049-16262). (a) Band 1. (b) Band 2. (c) Band 3. (d) Band 4. (e) Band 5. 

(0 Band 6. (g) Band 7. 


terize the scan -direction effect, the ratio of forward-scan signal 
level to reverse-scan signal level was computed as a function of 
scan angle (pixel number). We found that the droop had an 
exponential character and appeared to be proportional to the 
magnitude of the signal. This behavior is consistent with that 
of an ac coupled system with not quite complete dc restora- 
tion. To produce a preliminary model of the effect, it was 
convenient to fit the data to an exponential decay model. 
From these fits, we derived the following preliminary correc- 
tion model. 


1 + A [exp(-KM)~ 1 ] 

where 

S 0 (p) is the original signal for pixel “p”. 

S(p) is the corrected signal value for pixel “p”. 
p is the pixel number, counting from West edge of scene. 

Pf, p r are the pixel offsets. 

M = ( Pf + p) for forward scans = pixel number referenced to 
start of scan. 

(Pr ~ P ) for reverse scans. 


K is the time constant (reciprocal of the number of pixels of 
active scan required for 63 percent of the decay to occur). 

A is the factor determining magnitude of total decay. 

The model was fit with the constraint that A and K be iden- 
tical for forward and reverse scans. The result of fitting this 
model to data from Band TM1 from the Iowa scene is illus- 
trated in Fig. 12(a), and the predicted signal decay during 
forward and reverse scans is shown in Fig. 12(b). From this 
model, one can see that the decay begins prior to the first 
pixel, having decayed approximately 1 percent by that time. 
The total decay is approximately 2.25 percent by the end of 
the active scan. This results in the mean signal level of the for- 
ward scan being approximately 0.75 counts higher than the 
mean signal level of the reverse scan at the West edge of the 
scene, and approximately 0.75 counts lower at the East edge 
of the scene for Band TMl. 

For this case (Band TMl, Scene 40049-16262), the p f and 
p r required for the fit resulted in the exponential decay be- 
ginning at the time of dc restore (different times for the two 
scan directions). During dc restore, a dark (no radiance) level 
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is presented to the detectors, and a dc level capacitor is charged 
to a value which results in approximately two signal counts. 
This charge was designed to be held for the entire scan, having 
a time constant on the order of minutes [21]. The empiri- 
cally derived scan-direction effect model suggests that this 
charge may be decaying much faster (time constant of the 
order of 10 ms). 

We plan to continue our investigation of this effect, refine 
the preliminary correction model, and extend the analysis to 
other bands. 

C. Level-Shift Artifact and Recommended 
Correction Procedure 

1 ) Description of Artifact: When down-track profiles (mean 
signal versus scan-line number) for individual detectors (every 
16th line) from CCT-BT and CCT-AT data were examined, 
two patterns of level shifts were observed as one progressed 
down the frame. All detectors exhibited to some degree these 
step changes in mean signal across entire scan lines, although 
for most detectors the magnitude of the level shift was <0.5 
signal count. Also, it was found that these level shifts were 
correlated among detectors. Each detector exists in one of 
four noise states representing the combinations of the two 
forms or patterns of level shift and their respective phases. 

Form #1 of the level-shift artifact is exemplified by the 
first band, as illustrated in Fig. 13. TM1 Detectors 4, 12, 10, 
and 8 show the greatest effect, with level shifts of 2.2, 1 .8, 1 .0, 
and 0.75 signal counts, respectively, occurring at random and 
usually having several to many scans at a given level before 
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Fig. 10. Atmospheric scan-angle effects in MSS/TM Bands (Model 

results). 


shifting to the other level. A very noise-free definition of this 
pattern was obtained through analysis of nighttime frames of 
data from the reflective-band detectors. More than 6000 pixels 
were averaged to obtain each scan-line average value shown in 
Fig. 14. 

The pattern of the second form of level shifts is exemplified 
by Detector 7 of TM7, as shown in Fig. 15, again from night- 
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Fig. 13. Down-track profiles: Scan averages of scene data, before cor- 
rection, band TM1. 


Fig. 11. TM scan-direction effect, (a) Band 1. (b) Band 6. 
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Fig. 12. Scan-direction artifact and model for Band TM1. (a) Ratio of 
mean forward-scan signal to mean reverse-scan signal, (b) Predicted 
decay of signal with time. 
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Fig. 16. Scan averages (by detector) of calibration shutter data, band 

TM1. 


time data. This form #2 noise pattern is characterized by much 
more regular and frequent level shifts than form # 1 , often two 
mirror sweeps in each state. 

While we have observed the form #2 level-shift noise patterns 
in all of the frames we have analyzed, the form #1 pattern was 
absent in the most recent scene we examined (Grand Bahamas) 
and a portion (^ 3 ) of the earlier Augusta night scene (Fig. 14). 
Nevertheless, level-shift artifacts have occurred sufficiently 
often to warrant development of a correction procedure. For 
example, level shifts of the magnitudes observed in TM1 could 
seriously interfere with quantitative applications of TM data, 
such as in bathymetry [22] . Although we have not identified 
the specific noise sources in the sensor, we believe we have 
characterized the level-shift artifacts sufficiently to develop 
empirical methods for substantially reducing their effects in 
scene data. 

2) Recommended Correction Procedure: A correction pro- 
cedure requires a diagnostic measure to determine the level-shift 
noise state of each detector on each scan line, in addition to a 
method for applying the corrections. Several alternatives for a 
diagnostic measure were examined. The most promising uses 
data values obtained when the calibration shutter is blocking 
the input radiation and no calibration lamps are in the field-of- 
view. We obtained several CCT-ADDS tapes from NASA/GSFC. 
These tapes contain the desired samples of dark calibration 
shutter data in two data windows, each 24 or 28 pixels in 
width. We averaged the samples in the window that followed 
the dc-restore procedure for each detector for each scan line. 
The resulting down-track profiles again reproduced the level- 


shift pattern observed in scene data. Upon comparing the 
calibration shutter profiles for TM1 (Fig. 16) with the corre- 
sponding profiles from scene data (Fig. 13) and nighttime data 
(Fig. 14), one can see that while they exhibit more random 
noise than the nighttime data (^25 samples averaged versus 
>6000), they do not exhibit the scene-dependent fluctuations 
evident in the scene data. In addition to serving as a diagnostic 
of the system noise state, the dark calibration shutter values 
also give an estimate of the level shift magnitudes. 

Therefore, the correction procedure we recommend would 
compute the average of observations of the dark calibration 
shutter, after dc restore, and subtract the corresponding value 
from each pixel value on the scan line for the particular de- 
tector. Since this subtractive value will most likely be non- 
integer, incorporation of this offset would best be accomplished 
in the computation of the Radiometric Look-Up Table (RLUT), 
so as to eliminate an additional rounding/truncation step. This 
would, however, require an update of that table for each mir- 
ror sweep (every 16 scan lines). Even more frequent updates 
(e.g., for each pixel) would be desirable when correcting for 
the scan-angle effects discussed previously. However, this 
could impose severe implementation and operational problems. 

An example of the potential effectiveness of this correction 
procedure is illustrated by a comparison of the down-track 
profiles of Figs. 13 and 17. The scene content is apparent in 
the TM1 data of Fig. 13, as is the effect of low-frequency 
noise. Fig. 17 is the result of subtracting the shutter data (Fig. 
16) from the raw data (Fig. 13). Form #1 noise appears to be 
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Fig. 17. Scan averages of scene data, after correction. Band TM1 (ob- 
tained by subtraction of calibration shutter averages). 


eliminated, as have most differences among the down-track 
profiles of the detectors. In a similar comparison for TM7, 
form #2 noise was greatly reduced. 

It has been recommended to NASA that this procedure be 
seriously considered for incorporation into the standard radio- 
metric correction procedure for Landsat-4 TM data. Since 
the computation method used in this example was applied to 
the scan averages and therefore did not represent any added 
quantization effects that may result from a modified RLUT 
and pixel-by-pixel corrections, the recommended next step 
would be for NASA to analyze more recent data and, if 
warranted, implement a full correction procedure, test it on 
a trial basis on a number of scenes (including those we have 
analyzed) and carefully compare the new images and their 
statistics with those from the existing procedure. 

VII. Summary 

Engineering analyses of image data have been conducted for 
the Landsat-4 MSS and TM sensors. The image data were 
generally found to be of high quality and the TM gave evidence 
of several improvements over MSS, in its spatial and spectral 
characteristics. 

The MSS analyses provided good evidence that Landsat-4 
MSS data are of the same high quality as previous three 
Landsats, with comparable dynamic range and target responses, 
a good detector equalization procedure, and an accurate linear 
response to received radiance. Only two artifacts not noted 
in previous Landsats were found, and they were quantified. 
One is a geometric scan line length variation in Landsat “A” 
tapes, a variation that should not be a problem for the routine 


processing performed by the EROS Data Center to produce 
fully corrected “P” type data. The other artifact is a coherent 
noise effect in all bands, having an amplitude of less than one 
count. This noise causes a minor but distracting diagonal strip- 
ing pattern on images of relatively uniform areas and would 
require substantial processing to be routinely removed or re- 
duced. A well-defined linear relationship was determined to 
adjust Landsat-4 MSS signals so they follow a radiance re- 
sponse characteristic equivalent to that of Landsat-3. Based on 
less data, a simliar relationship was determined for Landsats-2 
and -4. 

Landsat-4 TM data were found to have very good spatial 
resolution and overall good radiometric data quality. Radio- 
metric equalization of detectors was found to be close to the 
theoretical limit of quantization error, except for two relatively 
low-amplitude artifacts. The first was a scan-direction artifact 
caused by signal droop during the active portions of each for- 
ward and reverse scan. This artifact was most pronounced in 
TM1 of the reflective bands, while the thermal band (TM6) 
had a unique pattern. The second artifact was a low-frequency 
scan-related level shift which occurred in two forms, leading to 
four system noise states. Again, this artifact was most evident 
in TM1. Initial correction procedures have been developed 
and recommended for further evaluation. The presence of 
scan-angle effects due to atmosphere characteristics and scene 
bidirectional-reflectance effects was noted and it is suggested 
that attention be directed toward development of normaliza- 
tion procedures. 

Upon comparing coincident MSS and TM data for one agri- 
cultural scene, high correlation values were found between MSS 
bands and the comparable TM bands. When spectral trans- 
forms were applied, Greenness variables for the two sensors 
were almost perfectly correlated. Their Brightness variables 
were less well correlated, likely due to contributions from TM 
bands covering regions not covered by MSS. Striking dif- 
ferences were observed in their ability to spatially resolve de- 
tail in this interesting agricultural scene. 
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Aiis i it \« i: Analyses o| the cliaraeteristies of Larnhat Thematic Mapper iTM) image data are 
described. and results are summarized Emphasis is placed on radiometric charaetcrizut ion. 
development of response models, and on comparisons between data from Eandsats 4 and 5. 
In general. the data quality was excellent; liim e\ er. some anomalies were found. 1 1 nee main 
topics are (a) systematic w itliin-sean-line signal droop use. <b> landont seancorrelated level 
shills, and (e) radiomelrie Ixigti.il amplitude) relationships between Lumbal - \ and l.andsat- 
5. rhe systematic clmop/rise effect was lonnd in data from both Eandsats 4 and 5. Daytime 
signals droop across the scan line while nighttime signals in the reflective bands rise across 
tlu* scan line*. The maiinitnde of tin* droop rise appears to be a function of tin* signal mag- 
nitude and average value of the signal throughout a scan cvcle Scan-correlated level-shift 
noise also was observed in data from both sensors, but with different patterns. Low-ampli- 
tude. low-frequency coherent noise effects also were measured The analysis ol simulta- 
neously acquired Landsat-4 and Landsat-5 TM data permitted a direct empirical comparison 
of the relative radiometric responses of their respective spectral bauds. Relationships be- 
tween their respective signal values were developed, and sensor dvnamic ramie consider- 
ations are discussed It was determined that multiplicative factors ramiimi from 0.9S7 to 
1.145 were required to convert tin* signal counts from Landsat-4 TM spectral bands to 
corresponding Landsat-5 equivalent signals Radiance v alues exhibited correspondimi dilfer- 
ences, pointing to residual errors in radiometric calibration. Low -lev cl clippimi was evident 
in the radiomcti icallv corrected Landsat-5 bands 5 and 7 data The temperature ramie cov- 
ered bv the full -S-bit data ramie of Tf PS-processed TM band 6 data vvas found to be* ap- 
proximately 2(K)°K to 340°K. not 200 k to 320 k as specified. 


INTRODUCTION 

S ince mi-: lai nc it of Landsat-4 in July 1982. nu- 
merous studies of the quality of Thematic 
Mapper (TM) image data have been performed 
under flu* auspices of the National Aeronautics and 
Space Administration s (NASA’s) Landsat Image 
Data Quality Analysis (LIDQA) Program. As part of 
this program, we have performed engineering anal- 
yses of Thematic Mapper image data with our efforts 
concentrated on radiometric characterization of the 
sensor. In general, we have found the data quality 
to he excellent. However, anomalies do exist in the 
data from both Landsat-4 and Landsat-5 TM. The 
analyses of Landsat-4 TM image data were previ- 
ously described in detail bv the authors (Metzler 
and Malila, 1983b; Malila et ai, 1984; Metzler and 
Malila, 1983a) and are summarized below. This 
paper concentrates on recent analyses of Landsat-5 


* This research was sponsored by the l T .S. National 
Aeronautics and Space Administration, Goddard Space 
Flight Center, Greenbelt, MD. under Contract N ASS- 
27346. 


TM data and comparisons of tlu* radiometrv of the 
two sensors. Specific topics covered arc (a) within- 
line droop, a phenomenon whereby the signal levels 
of the sensor change sv stemat icallv during tlu* active 
scan; d>) scan-con elated level shills, an ellect which 
raises or lowers the signal level of all pixels in a scan 
line* or set of scan linos; and (c) comparison of 
Landsat-4 and Landsat-5 radiometric corrections. 
Other analyses of TM data anomalies may he found 
elsewhere in this issue le g.. kiefler ct a /.. 1985). 

Within-Line Droop 

Earlier examination ol Landsat-4 TM average 
scan lines indicated significant differences in the sig- 
nals returned from the western edge of a scene com- 
pared to those observed at the eastern edge of the 
same scene. This effect was most apparent in the 
shortest wavelength spectral band (band 1), and was 
observed in all the spectral hands to some extent. 
A combination of bidirectional reflectance, atmo- 
spheric, and shadowing effects, as well as sun-view 
angle geometry can explain the effect observed. 
More careful examination of the average scan data, 
however, revealed a confounding ellect due to dif- 
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Fig. 1. Landsat-4 and Landsat-5 TM nighttime within-scan rise effect — band 1 


(emit sensor response characteristics related to the* 
direction of scan in hands 1 through 4. The scan- 
direction difference took the lorm ol a droop in 
signal with time during active scan, which appeared 
as a signal decrease with increasing pixel number 
for forward scans, and a signal increase with in- 
creasing pixel number for reverse scans. This effect 
was found in nighttime reflective-band data as well, 
but taking the form of a signal rise with time instead 
ol a droop. 


Scan-Correlated Level Shifts 

In Landsat-4 TM data, an effect was analyzed 
which changed the signal of all samples within a 
scan-line or group of scan-lines by up to 2.0 video 
quantum levels (l)\). The changes were aperiodic, 
occurring at random intervals with the level shifting 
during mirror turn-around time. All allected detec- 
tors shifted levels at the same time, with the level 
shifts following one ol two patterns (most detectors 
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Fig. 2. Example nighttime reverse scan signal levels — band 1 


exhibited both patterns, but one was dominant). 
Oik* pattern was exemplified by band I detector 4 
with a peak-to-peak amplitude of 2.0 l)N, the other 
by band 7 detector 7. These two patterns were la- 
beled form #1 and form #2, respectively (later 
labeled ‘tvpe l-l and * t \ pc* 4-7, respectix ely. In 
Barker (1984)). 

METHODS 

All analyses to characterize the radionietry ol the 
sensors were performed on digital computer com- 
patible tape (CCT) data. Several type’s ol (Cl data 
were used, representing various stages ol mound 
processing as well as calibration data. The analyses 
described in this paper generally were performed 
nil full-frame TM image data, both to characterize* 
full-frame effects and to take advantage ol the large 
data volume (approximately 37 million pixels per 


band per frame) to improve tin* <|iialilv ol tin* sta- 
tistics generated. 

Two primarv methods wore used to average the 
full-frame image data. In one ease, to examine scan- 
angle effects, average scan lines wore computed by 
averaging the columns ol pixels down the entire 
frame, 'lb analv/e scan-direction effects, these av- 
erage scan lines were stratified by scan direction, 
with the forward and reverse scan data being treated 
separately. The other tvpe ol averaging involved 
computing down-track profiles by averaging the 
rows of pixels across the entire Inline, thereby com- 
puting an average signal value* for each scan line. 
Each of these analvscs was performed separately for 
each band and eac h detector ol tlu* sensors. 

Earlier investigations In the authors (Malila ct 
dl.. 1984) demonstrated the value of using reflec- 
tive-band (bands 1-5 and 7) TM data collected at 
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night for analysis of sensor data anomalies. The 
sensor sensitivities are such that no scene radiance 
is recorded, so any variations in the data are due to 
sensor noise effects. We again made extensive use 
of nighttime data in the analyses described herein, 
processing the data using the techniques described 
above. 

Two techniques were employed to compare the 
radiometrv of the Thematic Mappers on Landsat-4 
and on Landsat-5, using a special data set which was 
collected simultaneously with both sensors. The 
first technique involved selecting a subimage of 
3. 504, ()()() pixels (1980 lines by 1800 pixels) from the 
Landsat-5 image and spatially registering it to the 


Landsat- \ TNI image data. This registration was per- 
formed to subpixel accuracv using 50 control points 
and nearest neighbor resampling. The data in each 
subimage* were averaged using 5 x 5 pixel cells to 
reduce any misregistration effects and to reduce the 
data volume while still retaining the* data diversity. 
Linear regressions were performed with data from 
the two averaged images. Multiplicative and addi- 
tive factors were computed for each band which can 
be used to relate the signals from one sensor to 
those expected from the other sensor for the same 
input radiance. 

For the second comparison, histograms were 
computed for the subimages described above, and 
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Scene 40037-02243 



Fig. 4a. Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 000 


Band 1 Forward Scans 
Scene 40037-02243 



Pixel Number 


Fig. 4c. Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 010. 


a histogram matching technique was employed to 
compute multiplicative and additive coefficients lor 
relating data from one sensor to that ol the other. 
Unlike the histogram matching procedure used iu 
TM ground processing which equalizes means and 


Band 1 Forward Scans 
Scene 40037-02243 



Fig 4b Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 001. 


Band 1 Forward Scans 
Scene 40037-02243 



Fig. 4d Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 011. 


standard deviations (Barker ct ai. 1983), a proce- 
dure Based on matching the cumulative distribution 
(unctions of the two data sets was employed. This 
technique is believed to have some advantages in 
making better use of the full range and distribution 
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Band 1 Forward Scans 
Scene 40057-02243 


Band 1 Forward Scans 
Scene 40037-02243 
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ig. 4e. Nighttime forward scan signal rise for scans 
receeded by calibration lamp state 100. 


2000 3000 4000 

Pixel Number 

Fig. 4f. Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 101. 


Band 1 Forward Scans 
Scene 40037-02243 



Fig. 4g. Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 110. 


Band 1 Forward Scans 
Scene 40037-02243 



0 1000 2000 3000 4000 5000 6000 


Pixel Number 

Fig. 4h Nighttime forward scan signal rise for scans 
preceeded by calibration lamp state 111. 


ol data values in the histogram. Additionally, the 
histogram matching approach has much less strin- 
gent registration requirements than pixel or region 
matching approaches. For a region with a very small 


perimeter/area ratio, effects of slight misregistration 
would he minimal. The cumulative distribution 
function gives percentage of observations having 
signal values less than or equal to the designated 
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Band 1 Forward Scans 



Increase in Scene Mean (DN) 


Fig. 5. Relationship of magnitude of signal rise and dif- 
ference between scan-cycle mean and scene mean. 


signal value*. Interpolations were made to obtain the 
signal values corresponding to integer percentage 
values. Excluding end points, a regression of the 
corresponding percentile signal levels from the two 
sensors provided the desired correction coefficients. 

RESULTS 

Within-Line Droop 

The single nighttime Landsat-5 TM scene (ID 5- 
0052-02182, Harrisburg. PA) available to us was 
used to quantify the vvithin-scan droop/rise ellect in 
Lmdsat-5 TM data. The average nighttime scan 
lines for band 4 of Landsat-5 TM are illustrated in 
Figure 1. along with data from Landsat-4 TM lor 
comparison, both forward and reverse scans are 
diown. The //-axes all have the same scale, i.e.. 0. 1 
l)\ full scale, to facilitate comparison between sen- 
sors. Note that for reverse scans, pixel position fiCHM) 
is sampled prior to pixel position I. I hcrelore, the 
effect is seen to be a signal rise* with increasing time 
lor both forward and reverse scans. In general, the 
within-scan rise has the same magnitude and time 
constant for the same band in each sensor. Magni- 
tudes are greater in daytime data and tlu* signals 
droop with increasing time, as will be discussed 
later. 

Baud 1 displays the greatest effect, with the mean 
reverse-scan signal increasing approximately 0. I 
DN during tlu* activ e scan. A simplified exponential 
decay model was fitted to the data lor each of bands 
1 through 4. For these bands, the time constant 


(time lor magnitude ol elfeet to deeav to 1 e of 
original value) which produced tlu* best lit ranged 
from 900 to 1 100 pixel sample times, lapproxi- 
matelv 0-10 milliseconds) for both Landsat-4 and 
Landsat-5 TM . 

The mathematical model used is expressed by the 
('((nation: 

St/0 = S„ + lie 1,1 (I) 

Where: 

S (/>) = signal returned by sensor for pixel p 
S„ = signal for p equal to infinity 
li ~ magnitude of total droop/rise 
7* = time (pixels) required for signal to 
change by (rl c /c of /i 

/ ) = pixel number, w ith count starting with 
lirst image pixel (west-most for forward 
scans, east-most lor reverse daytime 
scans). 

Since the magnitude and time constant of the 
nighttime w ithiii-scan rise are essentially identical 
for Laiids.it- I and Laudsat-5. we would expect the 
dav time droop elici ts to be similar also. During day- 
time data acquisition when signal levels were much 
higher, we observed in Landsat-4 TM data a cor- 
responding increase in the magnitude of the droop 
elicit. In a dav time band I scene ( 1 1) 4-0040- Ib2()2) 
which bad a scene mean of til. 0 DN. the magnitude 
of the droop was observed to be approximately 
minus 1.5 DN. with a time constant equivalent to 
approximately 000 pixels. At night, tlu* magnitude 
of the rise was <0. 15 D\. still with a time constant 
of 000 pixels for band 1 The mean scene level at 
night was 2.3 DN. Although qualitativclv the day- 
time I .andsat-5 effect appears similar to the dav time 
Landsat-4 effect, quantification of ibis effect in day- 
time Landsat-5 I M data awaits analysis of an appro- 
priate scene in which variations in seem* radiance 
have a rclativclv uniform spatial distribution. 

Hie droop, rise ellect was analv/ed further to es- 
tablish a hypothesis for its cause and a model Ibr its 
description and potential use in correction. While 
the magnitude of the effect does not appear to be 
strictly proportional to the scene mean, it does ap- 
pear as if the droop or rise is a drift toward the scan- 
cycle mean signal of the scene which also includes 
the signal values produced during shutter obscura- 
tion. calibration pulse, and IX! restoration. This 
scan-cv cle mean would be lower than the scene 
mean during the dav time due to the addition of the 
data acquired during shutter obscuration, and 
would be greater than the scene mean during night- 
time data acquisition, where the seem* itself is ef- 
fectively a continuation of the shutter obscuration, 
and (lie calibration pulses drive the scan-cycle mean 
to a level slightlv higher than the scene mean (see 
f igure 2). The hypothesis is that a-c coupling exists 
between the detector output and the analog-to-dig- 
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ital converter. producing a signal decay proportional 
to the departure from the scan-cycle mean. 

To test this hypothesis, a scene of nighttime 
Landsat-4 TM data (scene 4-0037-02243) was seg- 


mented based on the calibration lamp state ob- 
served by the sensors prior to each scan. As illus- 
trated in Fignre 3, the internal calibration lamps 
se(|uence through the eight possible states, re- 
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Fig. 6 Level shifts for Landsat-5 TM band 3 nighttime data 


main ini' in each state for approximately 40 scans (20 
forward/ re\ erse scan cycles). During nighttime re- 
flective- hand data collection, the signal pulses re- 
sulting from viewing these calibration lamps at the 
end of each scan are the only signals available to 
shift the scan-cycle mean from the scene mean. All 
scans with lamp state 000 (no lamps on) were 
grouped into one subimage, and seven other sub- 
images were created for the other seven lamp states 
(001, 010, Oil, 100, 101, 110, and 111. where each 
binary digit represents the state of one of the three 
calibration lamps). Average scan-lines were com- 
puted for each of these subimages, then smoothed 
and displayed as plots of mean signal level versus 
pixel position (see Figures 4a — 4h). Qualitatively one 


can see that the ( fleet is greatest when the calibra- 
tion pulse adds the most to the scan-cycle mean 
(state 111. all lamps on), and is nonexistent in the 
case of no calibration pulse (scan-cycle mean equal 
to scene mean). 

Quantitative support for the hypothesis of drift 
toward a scan-cycle mean was derived from data on 
the calibration tape (CCT-AD1XS) associated with 
the image data. From the CCT-ADDS. the magni- 
tude of the calibration pulse for each scan line could 
be computed, which in turn allowed calculation of 
the scan-cycle mean for each scan. The rise for each 
scan was computed and plotted against the differ- 
ence between scan-cycle mean and scene mean as 
illustrated in Figure 5. Regression analysis indicated 
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Table 2. Magnitude and Phase of Level-Shift Noise in Lanosat-5 Scene 5-0052-02182 
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3 

10 

0.257 

1 1.8 

< 

10 

-0 1 10 

8.0 


* Negative amplitudes indicate level shifts with phase shifts of 180 relative to band 3 detectors 


an excellent fit (R 2 of 0.99), which strongly supports 
the hypothesis, indicating that the parameter B in 
the model expressed above is a function of the dif- 
ference between the scan-cycle mean and the scene 
mean. Although this analysis was performed only 
for forward scans of band 1 of one Landsat-4 TM 
image, experience to date indicates that the result 


may be extended to both scan directions of bands 1 
through 4 of both Landsat-4 and Landsat-5 The- 
matic Mappers with a high degree of confidence. 

This droop/ rise effect has been observed for the 
primary focal plane bands only. For both Landsats, 
bands 5 and 7 show essentially no change in mean 
signal level within the scan line, with perhaps a 
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1980 by 1800 Pixel Subimage of Scene 4-0608-15463 
Landsat— 4 TM Band 7 
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1980 by 1800 Pixel Subimage of Scene 5-0014-15460 
Landsat-5 TM Band 7 
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Fig. 7. Landsat-4 and Landsat-5 TM band 7 histograms for coincident regions. 


slight change in tlu* opposite direction to that seen 
in hands 1 through 4. Band (i mean signal levels 
have been observed to change within scan lines in 
a variety of patterns. Detailed analysis of potential 


within-scan effects in band 6 is made more difficult 
b\ the absence ol any constant scene data compa- 
rable to the nighttime data in the reflective bands. 
Even a completely uniform ground scene would 
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1980 by 1800 Pixel Subimage of Scene 4—0608 — 15463 
Landsat — 4 TM Band 5 
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1980 by 1800 Pixel Subimage of Scene 5—0014 — 15460 
Landsat-5 TM Band 5 



have varying atmospheric effects in different parts 
of the scene. 

This droop effect should not cause serious prob- 
lems for most users. However, it can confound at- 
tempts to extend signatures from one side of a scene 


to the other and can introduce handing (stripes 16 
lines wide, or 17 lines in geometrically corrected 
CX.T-PT data) at the scene edges. Implementation 
of the proposed exponential model would require 
pixel-by-pixel correction and could prove costly in 
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terms of computation time. It is our understanding 
that NASA and NOAA will leave it to the individual 
users to determine the importance of correcting for 
this efleet and actually performing Jtlie correction. 

Also apparent in Figures 4a — lli are oscillations 
superimposed on the rise cllrct. These oscillations 
are coherent noise found in all rcllcctivc hands ol 
both Landsat-4 and Landsat-5. Although quite ob- 
vious in these plots derived from nighttime data, 
the peak-to-peak amplitudes are quite small (<0.75 
DN in unfiltered data, <0.05 ON in these smoothed 
plots) and have not been observed in daytime data. 
The cause of this approximately 400 Hz (202-204 
pixel period) noise is undetermined. 

Scan-Correlated Level Shifts 

In the Landsat-4 TM data we have examined, 
type 4-7 scan-correlated level shills are always 
present, and the signals often shift states with a reg- 
ular period. Sean-correlated shifts of type 4-1 are 
present in most, but not all data, and the typo 4-1 
pattern tends to remain in one state or the other for 
several scans of the scan mirror. I he peak-to-peak 
amplitude for each affected detector for each form 
of the shift is essentially constant in all cases where 
that form of the noise exists. The phase relationships 
between the affected detectors also remain constant 
in all images (i.c., band 7 detector 7 is alwavs in its 
high type 4-7 state w hen band 5 detector 8 is in its 
low type 4-7 state). Figure 14 of Malila rt ill. (1984) 
illustrated both patterns of level shift for the lb band 
1 detectors of Landsat-4 TM for a night scene. Kel- 
ative magnitudes and phases are readily apparent 
from the illustrations. Table l provides the quanti- 
tative results, giving the magnitude and phase of 
each level-shift pattern for the 90 Landsat-4 TM re- 
flective-band detectors. 

Initial analyses of Landsat-5 TM data indicated a 
similar effect, but with only one pattern (Malila and 
Metzlcr, 1984; barker 1984). We examined night- 
time reflective-band data to provide quantification 
of the magnitude and phase relationships of the ef- 


fect. f igure 0 illustrates the level shifts for band 3 
of Landsat-5 TM. tlu* band most affected by this 
noise. The plots were produced by computing the 
mean signal level for each scan for each detector of 
each band and by plotting these scan-line means 
versus the scan number. In these plots, the max- 
imum peak-to-peak amplitude is approximately 0.5 
DN. Table 2 contains the quantitative results for the 
re I lectiv e-band detectors of the Landsat-5 TM. It 
can be seen that nearly all detectors are affected, 
although the magnitude is very low (<(). 1 DN) for 
main. Band 3 shows the greatest effect although 
band 2 detector 1 is the single most affected de- 
tector with a level shift >0.5 DN. The compares 
with the maximum shift of 2.0 DN measured for 
Landsat-4 band l detector 4. Several detectors did 
not display any measurable efleet in this scene. 
They are: baud I detectors I, 3, 5. 9. 13, and 15; 
band 2 detector 1; band 4 detectors 8, 10, 12. and 
10; band 5 detectors 2. 4, 7. 10, and 13; and band 
7 detectors 1. 2. 5. and 15. As seen in Landsat-4 
TM data, patterns of phase and magnitude of the 
level-shift effect within a band often place the de- 
tectors into odd-even groups. As with the within- 
scan droop, the confounding effect of scene data 
prevents analysis of this type for band 6. For this 
hand, shutter data may be used to quantify any level 
shifts but with somewhat lowered precision. 

Although these level-shifts are strikingly ev ident 
in the nighttime reflective data, where the scene 
makes no contribution to tlu* observed signal level, 
they are of the same magnitude in daytime data and 
even there can cause noticeable striping. The mag- 
nitude of these shifts and the large number of scenes 
in which they occur places a high value on tlu* cor- 
rection of the effect for some applications. Fortu- 
nately. the constancy of the magnitude permits rel- 
ative! v simple correction techniques. Since the level 
shift remains constant for the entire scan, the shifts 
are also observable in the shutter data collected at 
the ends of each scan line. Based on this, several 
methods of correcting for level shifts have been pro- 
posed which appear effective in reducing the effect 


Table 3. Coefficients for Converting Landsat-4 DN to Landsat-5 DN (Scenes 4-0608-15463 and 5-0014-15460. 

15 March 1984) 


Band 

A 

Laiulsat-5 

li 

TM = A*(I 
SK. 

auidsat-4 TM) + li 
IF 

Bauge of Data Values 
(DN) 

L.mdsat'4 Lamlsat-5 

1 

1.0438 

-3.538 

0.131 

0.99943 

73-109 

73-111 

2 

1.1200 

-2.719 

0 134 

0.99922 

20-52 

26-56 

3 

0.9809 

-3.078 

0.142 

0.99975 

20-77 

22-72 

4 

1.0030 

-4.027 

0.078 

0.99995 

11-92 

7-88 

5 

1.1452 

-7.330 

0.100 

0.99999 

6-154 

0-109 

6 

1.0040 

-0.711 

0.119 

0.99956 

114-148 

113-148 

7 

1.0923 

-0.244 

0.054 

0 99999 

3-80 

0-88 


Note: If the DN computed for Landsat-5 is <0. sudstitute 0 If it is >255. substitute 255 
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Table 4. Coefficients for Converting Landsat-5 DN to Lanosat-4 DN (Scenes 4-0608-15463 and 5-0014-15460. 

15 March 1984) 


Hand 

A 

Landsat-4 

li 

TM = V 1 1 . 
s i; 

indsat-5 TM) + H 

IF 

Hangc of 1 )uta Values 
(ON) 

Landsat- 1 Landsat-5 

1 

0.9580 

4.490 

0 1 45 

0.999 44 

74 109 

74- 1 1 1 

2 

0.8928 

2 )27 

0 120 

0 99922 

20 52 

20 50 

4 

1.0142 

4.720 

0 114 

0 99975 

20 77 

22-72 

4 

0.9970 

4.01 1 

0 (478 

0 99995 

1 1 -92 

7 -88 

5 

0.8732 

0 401 

0.094 

0 99999 

0-15 4 

0-109 

0 

0.9900 

0.71 4 

0 1 fS 

0.99950 

III! IS 

114-1 48 

7 

0.9155 

5.717 

0.0 49 

0 99999 

4-80 

OSS 


Note If the DN computed for Landsat-4 is 255 substitute 255 


( Barker, 1984: Fischel, 1984: kognt ct at.. 1984. 
Malila ct at .. 1984; Mctzler and Malila. 1984a; 
Murphy ct at . , 1984). TIk* general approach is to 
detect the presence of the shift (normally l>\ looking 
at shutter data), then to subtract (or add) the know u 
magnitude of the shift to each pixel in the allot ted 
scan line. 

TM Landsat-4 Versus Landsat-5 
Radiometric Comparison 

Radiometric matching of the Landsat-4 and 
Landsat-5 TM sensors was lacilitated l>\ the a\ ail- 
ability ol a unique set of radiometricalb corrected 
data collected simultaneously by the two sensors 
and registered to subpixel accuracy as described 
above. Same-band images from the two sensors 
were very similar in appearance, although exami- 
nation on an image display system required dif- 
ferent gain and offset factors to be applied to achiex e 
identical brightness and contrast for each pair of im- 
ages. Since both the Landsat-4 and Landsat-5 
scenes were processed through Ties (Thematic 
mapper Image Processing System), it was expected 
that radiometrically corrected products would have 
essentially identical corrected signal values for the 
same scene viewed at the same time. In addition to 
multiplicative 1 and additive differences, clipping of 
the Landsat-5 data values was obvious in both bands 
5 and 7 at the low radiance end of the dynamic 
range. The band 7 low-level clipping is apparent 
from a histogram of signal-level frequency for band 
7 for both Landsat-4 and Landsat-5 (see Figure 7). 
The pixels with values zero to six in the Landsat-4 
scene are all mapped to value zero in the Landsat- 
5 seem*. Although the offset was nearly as large lor 
band 5 (see Figure 8), fewer data values actually 
were clipped (0.3 percent of the scene versus 4.2 
percent in band 7). 

As noted earlier, band-bv -band comparisons were 
carried out using two different techniques: (l) 
regression of signal values from the coincident pixels 
or regions, and (2) regression ol signal values asso- 
ciated with specific histogram percentile classes. 


When clipping was not present, either technique 
produced essentially the same results. Where clip- 
ping was present, regression of matched areas led 
to smaller additive terms and larger multiplicative 
terms, a result deemed erroneous after inspecting 
the histograms. For this reason, coefficients from 
the histogram matching approach are presented 
here, fable 4 presents the multiplicative and addi- 
tive coefficients to convert Landsat-4 TM signal 
levels to Landsat-5 TM equivalent values: and Table 
4 contains the coefficients to convert Landsat-5 sig- 
nals to Landsat- 1 values. It should be noted that 
while this was a simultaneously collected data set, 
and therefore nearly ideal for this type of analysis, 
the correction coefficients presented are valid only 
if ground processing parameters are not changed. It 
should also Ik* noted that data which have been 
clipped as in Landsat-5 bands 5 and 7 can not be 
retriev ed — all the /crocs in Landsat-5 band 7 data 
will be converted to sixes in Landsat-4 band 7, 
w hereas Landsat-4 band 7 would have recorded the 
same pixels with signal levels ranging from zero to 
six. In using these conversion equations, resultant 
DNs less than zero should be assigned the* value 
zero: DNs greater than 255 should be* assigned 255. 

Converting the pixel values to radiance levels via 


Table 5 Landsat-4 and Landsat-5 TM Radiance 
Conversion Parameters (Scenes 4-0608-15463 and 
5-0014-15460. 15 March 1984) 



Hadiuncc = 

.VO 4- W‘ 0\ (iiiW Kiir sr 

jim)) 



\() 

\l 



(ni\V true sr jim)) 

(mW/tcni" sr 

JJLIII ))/ 1 )\ 

Hand 

1 -mdsat • 1 

1 .andsat-5 

1 Mindsat- 4 

1 .audsat-5 

1 

-0. 15(H) 

-0 15(H) 

0.00021 

0.00024 

2 

-0.2802 

-0.2805 

0 I 1750 

0.11750 

4 

-0.1204 

-0.1 191 

0.08001 

0.08059 

4 

-0. 150 4 

-0 15(H) 

0.08145 

O.OSI44 

5 

-0.0472 

-0.0470 

0.010SI 

0.01081 

0 

0.1252 

0.1 248 

0.00509 

0. 00503 

7 

-0. 15(H) 

-0 15(H) 

0.00570 

O.(H)508 
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Table 6. Landsats-4 and Landsat-5 TM Regressions of Radiance Values (Scenes 4-0608-15463 and 5-0014-15460. 

1 5 March 1 984) 


Hattd 

\ 

1 .andsat- 

li 

5 1 M 
S 1. 

\ 1 1 .andsat - 1 TM ) ' 
R J 

li 

Range <il Radi 
0nV\ icnr 

1 atndsat • 1 

.nice Values 
sr (jliiiu 

Landsat 5 

1 

1 .0 13.5 

0 205 

u ouq 

0 999 13 

1 25 0 1 1 

1 23 (» 52 

2 

1 1 190 

- 0 2S5 

o u|() 

0 99922 

2.70-5. S7 

2.79 -0 2S 

3 

0.9-S03 

-0 297 

0 01 1 

0 99975 

1 .98 -ft. Oft 

1 07 - 5 OS 

4 

1 .0027 

-0 370 

OOOt) 

0 9999.5 

0.77-7.33 

0 12-0.9S 

5 

1 1 452 

-0 074 

o not 

0.99999 

- .03- 1 .03 

-.03-1.79 

ft 

0.9932 

-0 002 

0 001 

0 99930 

0.77-0.97 

0.70-0.90 

i 

1.0885 

-0002 

0 oof 

0.99999 

-.13-0.31 

- 15-0 3.5 


the coellicicnts provided in the Radiometric Cali- 
bration Ancillary Record of the Leader Kile associ- 
ated with each band of image data (NASA. 1 9S3 ) did 
not resolve the discrepancy observed between the 
two sensors. Table 5 lists the multiplicative and ad- 
ditive* coellicicnts extracted from tape headers and 
used in the conversion. Table 0 is similar to Table 3 
in that it defines conversion of Landsat-4 signals to 
Landsat-5 equivalent signals but in terms of radi- 
ance instead of signal counts. It is not known at this 
time why the radiometrically corrected data are not 
more closely matched. 

An additional discrepancy was noted between the 
previously published band 0 temperature sensitiv itv 
range and the range* implied bv the coefficients 
listed in Table 5. Csing these coefficients to convert 
the range 0.255 DN to radiance* gives a radiance 
range of* 0.125 to 1.575 mW/icnr sr //in), repre- 
senting an apparent temperature* range of* approxi- 
mately 200 to 340° K. not tlu* advertised 200 k to 
320° K. This causes an increase in the* temperature 
dillerence represented by a change of I 1)V 1 he 
specified 200°k t<>32()°k temperature range actually 
spans approximately 03-190 DN versus the speci- 
fied 0.255 DN. Kor Landsat-5 TM. the radiance 
range is very slightly diflerent (0.121 to 1.300 \\\W/ 
(enr sr //in)), still giv ing a range of apparent tem- 
perature of approximately 200 C K to 340 k tor a I)\ 
range of approximately 03-193 for apparent tem- 
peratures of 200° k to 320 k). Users unaware of 
these differences max incorrectly derive tempera- 
tures from TM band 0 data. 

SUMMARY 

Landsat-5 TM image data were found to be quite 
similar to Landsat-4 TM data, both in terms of high 
overall quality and in the presence of several anom- 
alies. Detailed analysis revealed a systematic 
vvithin-scan drift (or droop/rise) of the signal from 
the scene mean toward the overall scan-cycle mean 
in spectral bands 1 through 4. T he magnitude of this 
drift ranged from minus 1.5 DN (daytime) to +0. 15 
DN (nighttime), depending on scene content. The 
drift was fitted with a simple exponential decay 


model and found to have a time* constant equivalent 
to about one-sixth of a frame width. 

Scan-correlated level shifts are present in both 
Landsat-4 and Landsat-5 TM data. The maximum 
effect observ ed in Landsat-5 data was approximately 
0.5 l)N peak-to-peak. compared with a maximum 
of 2.0 DN observed in Landsat-4 data. I he level- 
shifts appear to be present in most if not all images, 
and effective correction procedures have been pro- 
posed. 

Although data from both Thematic Mappers are 
produced in radiometrically corrected form, com- 
parison of data acquired simultaneously by the two 
sensors revealed significant differences in their cal- 
ibration. In the reflective bands, the multiplicative 
factors required to convert Landsat-4 T NI data to 
Landsat-5 data ranged from 0.987 to 1.145, with 
corresponding additive terms of —2.7 DN to —0.2 
DN, and displayed evidence of low-level clipping 
in Landsat-5 bands 5 and 7. The thermal bands 
(band 0) were more closclx matched, but are cali- 
brated to have a lull-range temperature range of 
200 T k to 340° k instead of the advertised 200 k to 
320° k. 

REFERENCES 

Barker, f L. . 10-S I. Relative Radiometric Calibration of 
Landsat TM Reflective Rands: IASDSM-4 Science 
I nccsligalions Sinnnnm/. \ I. July. pp. 1 40- 1 SO. 
Barker. I L.. Abrams. R. B Ball. I) L.. and Leung, 
k C.. 1983. Radiometric Calibration and Processing 
Procedure for Reflective Bands on Landsat-4 Proto- 
flight Thematic Mapper: Proceeding v of Landsat-4 
Science Characterization Early liesults Sifinfiositnn . 
\ II. NASA Conference Publication 2355. p. 75 
Lischel. I) IBS I Validation of the Thematic Mapper Ra- 
diometric aifd (ieomelric Correction Algorithms: 
Il’.i'.i. Transact ions on Cleoseience anil lie mote 
Sensing . v. CL-22. no. 3. pp. 237-212. 
kietfer. II.. Cook. I). A . Lliason. K. M.. and Lliason. 
P. I . 19S5. Intrahand Radiometric Performance of 
the Landsat T hematic Mappers: Vhotogranonetric En- 
gineering anil licniotc Sensing, (this issue), 
kogut. J.. Larduinat. K.. and Fitzgerald. M . 198.3. An 
Vital) sis of New Techniques for Radiometric Coitcc- 


36 


ORIGINAL PAGE IS 
OF POOR QUALITY 


1330 NIOTOCKAMMICTMC ICNCIN 1C 1C HI NC 6c HICMOTIC SKNSINC. 1985 


tion of LANDSAT-4 Thematic Mapper Imam’s Ke- 
scarch and Data Systems. Inc.. Lanham. Ml) 

Malila, \V. A.. Metzler. M I) . Nice. D I’, and (.‘list, 
1C. I*., 1984. Characterization of I.ANDSAT-I MSS 
and TM Digital Image Data: IEEE Transactions on 
Geoscience and Remote Sensing. v. CIC-22. no. 3. pp. 
177-191. 

Metzler, M. I)., and Malila. W. A.. 1983a. Hadiometric 
Characterization of Thematic Mapper Fu II -frame Im- 
agery: in Proceedings SVSEJASP Conference on Tech- 
niques for Extraction of Information from Remotely 
Sensed Images, Hochester, NY. 

, 19831). Scan-angle and Detector K fleets in The- 
matic Mapper Hadiometry: Proceedings oj Eandsat-4 
Science Characterization Early Results Symposium . 
January 1983, v. II. NASA Conference Publication 
2355 pp. 421-441. 


Malila. \Y. A . and Metzler. M I) . 1984 Eighth Quar- 
terly Progress Report for Study of Sped rrd -Radio- 
metric (Characteristics of the Thematic Mapper far 
I .and Cse Applications ICnvironmciital Hesearch In- 
stitute of Michigan. Keport No. 1 H4000- 13- P. Oc- 
tolier. 

Murphy, J. M.. Butlin. T. Dull. P. F. , and 1*3 tzge raid. 

A. J.. 1984. Hevised Hadiometric Calihration Tedi- 
nicpie for LANDSAT-4 Thematic Mapper Data IEEE j 
Transactions on Geoscience and Remote Sensing, v 
CIC-22. no. 3. pp. 243-251. 

National Aeronautics and Space Administration. 1983. In- 
terface Control Document between the NASA Cod- 
darcl Space Flight Center (CSFC) and the Depart- 
ment of Interior FIU)S Data Center (1C DC) for 
Landsat-D Thematic Mapper Computer Compatible 
Tape (CCT-AT. CCT-PT). Hevision A: NASA. no. j 
LSD-ICD-105, |). fit). 


37 


N86- 166 


94 



INFRARED AND OPTICS DIVISION 


APPENDIX C 

Comparison of the Information 
Contents of Landsat TM 
and MSS Data 


PRECEDING PAGE BLANK NOT FILMED 


Comparison of the Information 
Contents of Landsat TM 
and MSS Data* 

William A. M alii a 

Environmental Research Institute ot Michigan. Ann Arbor. Ml 48107 


Aiisihu i A cninmiinicatinns-tlicnrv approach is taken to .uiah/e the dispersion and con- 
centration of signal values in various data spaces, irrespective of specific c lass membership. 
Entropy is used to (jtiantifv information, and mutual information is used to measure tin* 
information represented l>\ subsets of spectral variables. Several different comparisons of 
information content are made. These* include* comparisons of svstem design capacities, of 
data volumes occupied In agricultural data in the spaces defined In original bands and In 
transformed spectral ( lasseled ( apt variable’s, of the* information contents of original bands 
and Tasseled C lap variables, and of the information contents of I’M and MSS for the given 
agricultural data sets Also, the effects of sample* si/e. scene content, and (pianti/ation level 
are examined. 


INTRODUCTION 

I N ANALYSES OF Mt'LTISFFCTHAL cliltil Sets produced 
by imaging remote sensing systems, needs arise 
for comparing the amounts of information provided 
by individual spectral bands, by various combina- 
tions of bands, and by different sensors. Measures 
based on classification performance or signal vari- 
ance (e.g.. principal component analysis) are com- 
monly used for such comparisons. Classification pro- 
cedures require knowledge of the identity of the 
scene elements being imaged and usuallv involve 
assumptions on the form of the signal distributions 
and parametric descriptors of those distributions. A 
class-independent and nonparametric measure of 
information content can be described in informa- 
tion-theoretic terms and is used here to analyze and 
compare digital image data from the Landsat Mul- 
tispectral Scanner Subsystem (MSS) and from the 
Thematic Mapper (TM). 

C. Shannon (1948) developed entropy measures 
of the information content of communications sig- 
nals. Trice (1984) and Bernstein, rt al. (1984) made 
entropy calculations and comparisons of Lamlsat 
data on a band-by-band or component-by-compo- 
nent basis. Malila (1984) developed a procedure that 
takes into account dependencies among spectral 
bands and applied it to original and transformed ver- 
sions of Landsat data; those results are summarized 
herein and extended to include additional data sets 
and other considerations. 


* This research was sponsored by the l T .S National 
Aeronautics and Space Administration, (hxldard Space 
Flight Center. Creenbelt. Ml), under Contract NAS5- 
27346. 


METHOD 

A communications-theory approach is taken to an- 
alyze the dispersion and concentration of signal 
values in various data spaces. Entropy, as defined 
by Shannon (1948). is used to quantify information. 
The process of selecting a subset of bands is viewed 
as the transmission of data through a communication 
channel in which loss of information may occur, and 
the mutual information between input and output 
is used to measure information transfer, i.e., the 
information represented by the subset. 

Several different comparisons of information con- 
tent are made. These include (l) comparison of TM 
and MSS system-design information capacities, (2) 
comparisons of the TM and MSS data-space vol- 
umes spanned by the agricultural data in the spaces 
defined by both original bauds and transformed 
spectral (Tasseled Cap) variables, (3) comparison of 
the agricultural information content of original 
bands to that of transformed variables, and (4) com- 
parison of the agricultural information content of 
TM data to that of MSS. The effects of sample size 
and varied scene content are examined, as is the 
effect of coarser quantization. 

Basic Information Concepts 

Shannon defined self information, /(.v,), as a mea- 
sure of the information associated with knowing the 
occurrence of a signal state .v, which occurs with 
probability P (a*,); 

I(x,) = log, ~ -lofePCv,) (bits) (1) 

The rarer the event, the greater is the uncertainty 
about when it will occur and, consequently, the 
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greater is the information conveyed when it is ob- 
served. Entropy, given the symbol H. is the value 
of self information when averaged over all \ pos- 
sible states of x: 

v i 

H(x) = 2 P(\) log, — (2) 

» = l P(X t ) 

Entropy is at its maximum when all states or cells 
are equally likely. It can be reduced by decreasing 
the number of cells occupied, by hav ing a nonuni- 
form distribution or a concentration of observations 
in the occupied cells, or both. 

With two variables, the use of joint and condi- 
tional probabilities is necessary: 

H(xjj) = H(x) + H(y\x) (3) 

since 

P(xjj) = P(x)P(t/\x) (4) 

In computing the conditional entropy, the 
weighting assigned to each information term is the 
joint probability of the states involved, i.e.. 

N v N If . 

H(x\y) = S E P(x r y,) !<>gj — - — (5) 

i = i j = i Pl-v.l!/,) 

If we consider x to be the input to a communi- 
cation channel and ij to be the* output, we can define 
the mutual information transferred between them, 
i.e., / w (A\f/), as 

/ A/ (.V.(/) = //(A*) - //( A*| If) (C)) 


become multidimensional vectors X and 1. with X 
= X s ) and V = (V,,V 2 ,. . .A\L Usu- 

ally. X lf < \ t The information transfer achiev ed by 
the communication channel is used lien* in a general 
sense, to represent both simple selections of spec- 
tral band subsets and more complex transforma- 
tions, such as the Tasseled Cap transformation. 


Absolute vs. Relative Information Content 

Multispectral sensors produce signals that have a 
fixed maximum number of signal levels in each spec- 
tral band, usually expressed as a number of bits, 
e.g., six bits for f>4 levels in telemetered Landsat 
MSS bands and eight bits for 25b levels in Landsat 
TM bands. When the probabilities in the entropy 
equations are based on all possible combinations of 
those levels, absolute information measures will re- 
sult. These would, for instance, be appropriate 
w hen absolute radiometric calibration of data is uti- 
lized. 

Most current uses of multispectral data, however, 
employ techniques that utilize only relative ampli- 
tude information between signals from various 
scene elements. In these instances, the information 
resides in the number of spectral cells that are oc- 
cupied and the distribution of spectral values within 
them. Malila (1984) developed an expression that 
gives a relative entropy value. //,,. for any given 
data set. in terms of counts of occurrences of obser- 
vations in cells of the spectral space. It is repeated 
here (for six v ariables): 


//,<(*> = 




ilk hint 


Informa- 
tion if each Information loss due to con 
observation centration of the observa- 
vvere in a tions into a subset ol cells 
unique cell 


(Ta) 


This equation shows that the mutual information ex- 
changed is the difference between //(.v). the infor- 
mation content ol the input, and //(.\jf/), tlu* infor- 
mation loss or uncertainty about x when we are 
given the output //. When tlu* total information is 
transferred, //(.\j//) = 0 and l Sl (x:ij) = ll(x). At the 
other extreme, when // does not contain any infor- 
mation rclatable to .v. H(x\y) = //(.v) and therefore 
7 a/ (.v,i/) = 0, i.e., there is no mutual information. 
Figure 1 presents a concise graphical summary of 
these quantities and their interrelationships. Nu- 
merical examples are given by Malila (1984). 

Multispectral Extension 

The above concepts can be extended to multi- 
spectral situations by letting the variables x and (/ 


where: C’„ A/|||IJ is tlu* count of occurrences in the 

cell having Level i in X r Level j in A>, etc., 

and: is the total number ol observations in 

the data set being analyzed. 

More briefly, 

//«(*> = (71.) 

It also is informative to div ide the total informa- 
tion loss due to spectroradiometrie concentration of 
signals (from Equation 7) into two components, one 
due to the reduced number of spectral cells which 
are occupied (below the total possible) and the re- 
mainder which occurs w hen the duplicate observa- 
tions are not uniformly distributed among those 
cells, i.e.. 


41 


TM AM) MSS DATA IM'OHM \IIO\ COM KM 


1451 


JOINT ENTROPY 

H(x.y) 



ENTROPY 

H(y) OF 

OUTPUT 


CONDITIONAL l H(x}y) V i( X y) V H (yfr) CONDITIONAL 
ENTROPY (LOSS) ENTROPY 


Fig. 1. Summary of Information Relationships. 


H, 


- K. ii + I*, 


(•S' 


‘Kins *'Yrll _r I- 'mnl 

where: L K . r n is the cell loss or loss in number of 

cells, i.e., 

Cm = ll 'fe - I‘>k 2 .V„.ik = -l<»«2 ( ) 

V * ' nils ' 

and: L is the uniformity loss, 

L , = IL - 


Spectral Band Subsetting 

The selection of subsets of spectral bands is a spe- 
cial case of the mutual information expression, 

/ W (X;Y) = H(X) - H(X\Y) 

where V now is a subset. A", of the X variables, so 

/ A/ ( A\X') = //(X) - //(X|X") 

Whenever a variable, say X ;1 . is retained, its con- 
ditional probability term becomes unity, its contri- 
bution to /7(X|X') is reduced to zero, and its infor- 
mation content is retained as mutual information. 
Whenever a variable, say X , is eliminated, there is 
a loss of mutual information. This loss is represented 
by the conditional entropy term through all condi- 
tional probability components in which X (/ occurs on 
the left-hand side of the conditional probability in- 
dicator line but not on the right-hand (or given) 
side. 

Spectral Transforms 

Spectral transformations wore obtained by ap- 
plying the* linear-combination Tasselod (lap 
(iasca i*) transformations to MSS (kauth and 
Thomas, 197b) and six-band TM (Crist and Cicono. 
1984) data. The principal r\s< \r variables are 
Brightness and Greenness. The Brightness variables 
are positively weighted sums of all bands and re- 
spond to general changes in overall scene ref lec- 
tance. The Greenness variables are essentially con- 
trasts between near-infrared wavelengths (where 
healthy vegetation is more highly reflecting than 
soil) and visible wavelengths (where healthy v ege- 


tation tends to be less reflecting than many soils) 
and respond to the amount of vegetation present. 
These two v ariables capture 95 to 98 percent of the 
variability in MSS data from tv picul I’.S. agricul- 
tural scenes, while a third variable, called W etness, 
has been found to be significant in similar TM data 
(Crist and Cicone. 1984). The Tasseled Cap vari- 
ables. though related to principal-component vari- 
ables. have advantages over them in that the Tas- 
seled ( lap directions do not van with the scene con- 
tent. and they have more consistent intcrprctability. 

Also, principal-component analysis was utilized to 
obtain a different set of spectral variables for one 
comparison. All transformed values were rounded 
to the nearest integer before being analyzed. 

Quantization Effects 

To explore the influence of quantization on the 
resultant information content, the amplitude values 
were* requantized several times. At each step, the 
number of original digital counts per modified am- 
plitude* interval was doubled, thereby compressing 
the* data and reducing the number of bits per 
channel by one lor each step. 

DATA SET 

MSS and six-band TM data of two types were 
analyzed. T hese are (1) real Landsat-4 MSS and TM 
data acquired simultaneously from an agricultural 
scene in North Carolina and (2) data values synthe- 
sized from field-measured reflectance spectra of ag- 
ricultural crops and soils using an atmospheric 
model. T hese data were used in prior comparisons 
of the spatial and spectral characteristics of Landsat 
TM and MSS data(Malila, ctal. 1984; Crist, 1984). 
In the synthetic data, samples are primarily from 
vegetation at a variety of ground cover percentages, 
with many fewer examples of bare soil. All analyses 
of TM data are limited to the six ref lective bands; 
the thermal band was not analyzed in this effort due 
to its coarser spatial resolution, its dependence on 
emissive rather than reflective characteristics of 
scene materials, and lack of a comparable simulation 
data base. The TM frame was acquired on 24 Sep- 
tember 1982 and included a wide range of agri- 
cultural crop conditions, ranging from bare soil to 
green and senescent vegetation to crop residues. It 
also included some samples from water and vege- 
tation along the Atlantic Coast and from deciduous 
and coniferous trees. 


RESULTS 

Spectral Data Volumes 

Tin* diagram in Figure 2 helps describe the var- 
ious terms used here to designate spectral data- 
space characteristics, while Table 1 quantifies many 
of the observed values. Figure 3 presents informa- 
tion measures for several of those quantities, as a 
function of the number of data v ariables. First, the 
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Fig. 2. Illustration of Various Spectral Data Volumes. 


system-design capacities of the Laudsat-4 TM and 
MSS arc presented in terms of the number of bits 
transmitted to the ground and/or recorded on com- 
puter-compatible tapes (CCTs). For TM, the 
number of bits recorded on CCTs is the same as 
that transmitted (8 bits/channel). For MSS, how- 


ever, the six-bit telemetered data are expanded to 
seven bits on the CCTs, with only an apparent gain 
of information. Nevertheless, many comparisons in- 
volving MSS will use seven-bit data since that is the 
form in which we received them. For some others, 
a degradation to six bits was performed before anal- 
ysis. The greater information potential of the TM 
system design (reflective hands), us compared to the 
MSS system, is (juantilied as 48 vs. 24 bits in tele- 
metered data. 

Figure 3 also portrays the hvpercube volume or 
data-space volume spanned by the TM and MSS 
data of 'fable la. These volumes are computed by 
summing the bit equivalents of the observed data- 
value ranges (max — min + 1) in each band being 
considered. Upon comparing the fractions of their 
total data-space volumes that are spanned by data 
from the agricultural scene, one observes that the 
TM data fall nine bits short of capacity while the 
MSS data fall approximately six bits short of ca- 
pacity. 

Actual data dispersion volumes or relative entro- 


Table 1. Information Comparison for MSS and Six-Band TM Data Sets 

A. Values for Heal Agricultural Data 

(A common area from the \. Carolina scene) 


MSS 

Si\-KI;uid TM 

TM Gain 

iiiiiiiIkt hits 

tmmher hits 

hits 


Svstem capacity: Sensor 

0.17 x 10 s 

24 

0 2S x 10 n 

48 

24 

CUT 

0.27 x 1()‘‘ 

28 

0.2S x 10 n 

48 

20 

Ilvpcrcuho vol.. Sensor 

0.44 x 10“ 

18.7 

0. 13 x 10 1 - 

38.0 

20 

CUT 

0.32 x 10- 

21.6 

0.43 x 10 1 - 

38.0 

17 

Data dispersion pattern. 






• # Observations. ll„„ n 

3. IBS 

11.8 

13.015 

13.7 

1.91 (Spatial) 

• # Unique cells 

2.808 

— 

12.903 

— 


\!S^ fYT * Relative Entropy. ll H 

‘ , | • Entrnpv loss clue to 

/ hits per hand . . 

1 spectral concentration. 


11.4 


13.7 

2.27 (Total) 

Hi.,, 

— 

0.38 

— 

0.02 

0 30 (Spec ti ah 

MSS Sensor. • # Unique cells 

1,730 



12.903 




6 bits per band • Relative Entropy. II l{ 

— 

10.3 

— 

13.7 

3.34 flblal) 

• Entropy loss clue to 






spectral concentration, 
///<»>< 

— 

1.45 

— 

0.02 

1 . 43 (Spectral) 

B. Values for Synthetic Agricultural Data 






(Assumes equal spatial resolution) 
“System capacity (MSS: 0 hits/haud) 

017 x 10 s 

21 

0 28 x |0 n 

IS 

21 

Observed hvpercube volume 

0. 10 x 10 : 

20 

0.99 x 10‘- 

10 

20 

Data dispersion pattern: 






• # Observations; ll„ MX 

2.270 

11.15 

2.270 

11.5 


• # Unique cells 

817 

— 

2.2(jO 

— 


• Relative Entropy. ll H 

• Entropy loss clue to 

— 

8.94 

— 

11.14 ' 

I 2.20 

spectral concentration. 

i/ 

— 

2.21 

— 

0.01 , 

| (Spectral)* 


^ 

' (TM gam over seven-bit simulated MSS data was one bit ) 
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Fig. 3. Comparison of Landsat TM and MSS Informa- 
tion Capacities. 


pies (see Figure 2 and Table 1) were found to be 
substantially smaller than the hyprreubr volumes, 
owing to correlations between bands and the limited 
numbers of observations. Results for the real TM 
data are shown in Figure 4 and lor both TM and 
MSS (7 bits/band; CCT) in Figure 3. Note that these 
relative-entropy values lor actual information are 
substantially smaller than those reported by Price 
(1984) for similar comparisons in w hich the sum ol 



NUMBER OF VARIABLES 

Fig. 4. Thematic Mapper's Utilization of Data Space. 


band values was treated as the joint information con- 
tent. T he data dispersion volumes in Figure 4 are 
measured by the relative entropies ol the best vari- 
able combinations, and represent tlu* relative infor- 
mation present in those sets. Most ol the informa- 
tion is contained in the first two or three variables, 
both the best and worst combinations are shown lor 
each system in Figure 3. T he number ol observa- 
tions analyzed establishes a maximum limit on each 
relative entropy v alue. As shown earlier in Filiation 
(7), the concentration of multiple observations 
(pixels) into indiv idual spectral cells reduces the in- 
formation content below the potential maximum. 
Table I shows very little tendency for TM pixels to 
do this, owing to the very large system capacity, 
spectral diversity, and fine gradation ol the TM 
bands. T he* MSS data show definite tendencies for 
multiple observations in spectral cells. 

Table I shows that the TM data represent 3.3 bits 
more information than the MSS sensor data, with 
approximately two bits being associated with spatial 
resolution (pixel size and number) and the re- 
mainder with spectral bands and radiometric reso- 
lution. Since the synthetic data have the same 
number of observ ations for both TM and MSS. they 
can be considered to have ecpial spatial resolutions. 
Thus, the 2.2-bit difference must be solely due to 
their spectral and radiometric properties. 

The above results were for a systematic sample 
from a larger area. 900 lines by 1300 TM pixels in 
size (430 x 030 MSS pixels). To explore the effects 
ol sample size and scene content on the information 
measure, the area was divided into nine subareas, 
containing varied types and amounts of the scene 
classes. W hen all 1.17 million TM pixels were in- 
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Fig. 5. Range of Information in Subsets of Bands. 
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eluded in the analysis (Data Set TM-C), an infor- 
mation content equal to 18.4 hits of the possible 
20.2 hits was computed, as shown in Table 2. For 
the corresponding 0.29 million MSS pixels. 13. S hits 
of the possible 18.2 hits were present as informa- 
tion. The two hits difference between maximum po- 
tentials is due* to the greater number ol I’M pixels, 
deductions below the maxima are due to reduced 
numbers ol distinct spectral cells and nonuiiilormitx 
of the cell populations. Bit equivalents of those 
losses are also indicated in Table 2. It can be seen 
that substantially greater losses occur for MSS data 
than for TM data, leading to a total difference of 
4.6 bits between the two data sets. 

Values also were computed using all pixels in each 
subarea. Mean values are given in Table 2 (Data Sets 
TN1-B and MSS-B) along with standard deviations 
to indicate the amount of variability found among 
the different scene areas. On the average, both 
types of losses are reduced from those found in the 
total data set, but variability among subareas is sub- 
stantial. Even smaller subsets of data were obtained 
for analysis by taking every tenth pixel in each sub- 
area; the average's and standard deviations of those 
values also are listed in Table 2 (Data Sets TM-A 
and MSS-A). For these. the loss of information by 
TM is very minor (0.20 bit), but the losses for MSS 
remain greater (about one bit). \\ barton ( I OS P sim- 
ulated TM and MSS data sets, analv/ed histograms 
of various sample sizes, and computed ratios of dis- 
tinct to total number of samples. Comparable num- 
bers were computed from the* average cell loss 
v alues and are giv en in the last column of Table 2 
as the percentage of cells which are distinct. The 
percentage for the largest real MSS data set is quite 
comparable to that for the largest set examined by 
W harton, who found 27 percent distinct cells among 


230. 100 samples. II is 50 percent for 28,800 samples 
and 85 percent for 3.000 samples are higher than 
the respective 34 percent in fable 2 for 32.300 sam- 
ples and 00 percent for 3.300 samples. For TM, 
W barton found nearly 100 percent distinct cases for 
even the 230. 400-sample case, versus 01 percent 
here for 130.000 samples, but lie considered seven 
rather than the six dimensions analyzed here and 
inc luded onlv samples from nine scene classes. The 
TM data in Table 2 retain much more distinctness 
than MSS as the sample size is increased, with SO 
percent distinct for 13.000 samples. 01 percent for 

130.000 samples, and 38 percent distinct for 

1.170.000 samples. 

Spectral Transformations 

Figure 4 also compares the data-space volumes 
spanned by original bands and Tasseled Cap trans- 
formed versions of signals from the agricultural 
scene (Table I \ sample). 'Three fewer bits per pixel 
arc* required to provide* the same* information using 
the transformed variables than would be required 
by the* original hauels. I bis ellect potentially could 
be useel to reduce* telemetry requirements; differ- 
ences might be even greater for data sets with a 
broader range of sce*ne amplitude's. 

For the synthetic* MSS data s<*t. a comparison was 
made of the* information content ol original baud 
values and two tvpe*s e>l transformed variables. 
T \s( \ p variables and principal-component vari- 
ables. They were found to be essentially identical. 
The* equality of llu* complete sets of variable's is in 
keeping with theoretical considerations ol linear 
transformations. 

To compare* with the* original-band values of 
f igure 3. relative* entropy values lor the* best and 
worst l \s< \r subse ts e>f each size* are presented 


Table 2. Effects of Sample Size and Scene Diversity on Information Content 


Data 

Sot 

Number 

of 

Pixels 

II 

Maximum 
Possible 
Kelatix t» 
Entropy 

(bits) 

n„ 

Actual 
Kelatix e 
Kutrnpx 

tbits) 

/...« 

1 .oss in 

Number 
ol Cells 
(bits) 

k»ni/ 

luiformity 
1 .oss 
tbits) 

— — X 1(H) 

■V 

Pent *nt 
1 )istinet 
Cells* 

TM-A 

1.3 x 10‘ 

13.67 

13.41 * 

0.162* 

0.091* 

89. 1 




(0 21) 

(0 136) 

<0.078) 

|73-99| 

TM-B 

1.3 x 10 5 

16.99 

15.66* 

0.711* 

0.615* 

61.1 




1 1.37) 

(0.675) 

(0.704) 

|.J8-90| 

TM-C 

1.17 x 10“ 

20. 16 

18.41 

0 791 

0.93) 

37. 8 

MSS-A 

3.3 x 10 1 

11.69 

10.72* 

0.604* 

0.361* 

63.8 




(0.47) 

(0.331) 

(0 1 IS) 

[)4-S8| 

MSS-B 

3.23 x 10' 

14. 99 

12.27* 

1.339* 

I 179* 

34.4 




(10 4) 

(0.093) 

tO 380) 

1 16-611 

MSS -C 

2.93 x 10’ 

18.10 

13.81 

2.149 

2.2(H) 

22.5 


* Denotes mean of values from nine subareas 
( ) Denotes standard deviation of those values. 

+ Computable from average bits of cell loss, i e 100 * 2 exp( - l c#ll ) 
( ] Denotes range of values computed for individual samples 
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in Figure 6. In this ease, we find an even greater 
disparity between best and worst combinations, 
owing to the decreased information content of the 
last TASCAP variables. Here again. relati\el\ little 
information is gained by the inclusion of more than 
three variables. 

Dimensionality 

Figure 7 displays relative entropy values com- 
puted for the first three Tasseled Cap components 
of TM and MSS data from the agricultural scent* 
(Table l A sample). (The MSS data were in ( .’( I form 
at seven bits/band.) The first three components are 
individually quite similar for TM. but there is a sub- 
stantial decrease (3.3 bits below brightness) for the 
third component of MSS (Yellowness). This is con- 
sistent both w ith many investigators experiences in 
finding MSS data of agricultural areas to be pri- 
marily two dimensional and with recent studies 
which have found a substantial amount of iiilnrma- 
tion in the TM Tasseled Cap Third Component 
(Crist and Cicone, 1984). Throughout this compar- 
ison, TM values are greater than the corresponding 
MSS values, for example the TM brightness value 
is 6.7 bits compared to 5.8 bits for MSS. 

When pairs of components are considered, we see 
substantial increases in total information, as would 
he expected with the addition of a second variable; 
the value for TM Brightness/Greenness is 4.8 bits 
greater than for Brightness alone, and the corre- 
sponding increase for MSS is 3.7 bits. However, 
differences do appear between MSS and TM. 
Whereas the value of the brightness/Creenness pair 
for MSS is substantially greater than the other two 
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Fig. 6. Range of Information in Subsets of Tasseled 
Cap Variables. 


(approximately two bits greater than Creciincss/ 
Third Component), there is relatively little differ- 
ence* (less than 0.4 bits) among the three pairings 
from TM data, pointing to a higher dimcnsiniiuiilv 
in TM. 

T hree components captured the vast majority of 
information for both systems. I low ev or, tlu* fact that 
the gain in going from two to three components was 
nearly as large for MSS (1.25 bits) as for TM (1,70 
bits) was somewhat surprising in view of the pre- 
v iouslv discussed two-dimensional character of MSS 
data. Furthermore, principal-component analvsis of 
MSS data showed nearly total representation of vari- 
ance by tlu* first two components. The MSS gain 
likely is due to the brightncss/Grceimcss plane 
hav ing a thickness of several counts in the third di- 
rection. even though this third component was un- 
correlated w ith the others. T lie observed values also 
indicate (hat differences do exist among these var- 
ious measures ol multispcctral signal properties. 
The TM data pattern also may be somewhat planar 
in three space, although not aligned as well with 
any component axis; correlations with the Third 
Component were —0.09 for brightness and 0.30 for 
Greenness in this data set. None of these observa- 
tions. however, should diminish the utility of Tas- 
seled Cap transforms for physical interpretation of 
data values and agricultural scene characteristics. 

Noise 

Noise in multispectral data was not considered 
explicitly in the results presented thus far. Sensor 
noise effects certainly were present in the real 
Landsat data, and natural variations of crop obser- 
vations were present in both the real and sy nthetic- 
data. Noise can add variance to signals and increase 
the number of spectral cells occupied (above that 
for no noise), thereby creating an apparent infor- 
mation content greater than the true information 
content of ideal, noiseless signals. To explore suc h 
ellcets. the number ol discrete levels present in the* 
data sets was reduced by apply ing several different 
((iianti/ation factors (greater than unity ) to each 
band and computing the reduced information con- 
tent. The results are summarized in Table 3 for 
three subareas which had (relatively) high, medium, 
and low information content, respectively. The TM 
still had more information when degraded to seven 
bits per band but, by the time the amplitude data 
were degraded to five* bits per band, there was little 
difference between the corresponding TM and MSS 
data sets. 

SUMMARY 

An information-theoretic measure was defined 
and used to compare Landsat MSS and TM multi- 
spectral data. The measure quantifies signal disper- 
sion patterns, independently of class membership 
and distributional assumptions. It provides an alter- 
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Fig. 7. Comparison of Information Contents of TM and MSS Tasseled Cap Variables. 


is 


naU* method (to classification) ol measuring the ex- 
tent to which subsets ol bands or transformed vari- 
ables represent the total pattern. The relative en- 
tropy value is limited by the number of Obsei v ations 
being analyzed. Since* results do vary with scene 
content, analysts should insure that data sets being 
analyzed are representative of the problems under 
consideration. 

A number of observations were made. The TM 
system-design information capacity is much greater 
than that of MSS. The potential information capac- 
ities and the signal hypercube volumes of agricul- 


tural data were much larger than the information 
actually represented by signal dispersion patterns 
in the sets of data values analyzed. Tasseled Cap 
transformations preserved the information in orig- 
inal bands and offered a modest savings in bits over 
those original bands, a lad which might be iiselul 
in data compression approaches. Relatively lew 
multiple occurrences of spectral observations were 
found in the TM data sets compared to MSS. an- 
other indication of TM s liner partitioning of spectral 
space. For the* best' combinations of variables, rel- 
ative entropy magnitudes were more a function of 


Table 3. Effects of Quantization Detail on Infommation Content 


Sensor 

Number 

of 

Pixels 

Helativ e 
Scene 
( amiplexilv 


Relative Kntropv ihitM 
lor Indicated Number ol lilts per 

Hand 


S 

7 

0 

3 

1 

3 

TM 

1.3 x Id' 

High 

10 0 

13 o 

12 3 

0.0 

0.2 

4.0 


1.3 x Id'' 

Medium 

10. 1 

13 1 

in n 

7 .-> 

1 *S 

3.4 


1.3 x Id' 

Low 

13.3 

1 IS 

S3 

3.S 

3.7 

2.7 

MSS 

3.3 x Id 1 

High 

— 

13 S 

in o 

0.1 

0 . 1 

4.3 


3.3 x Id 1 

Medium 

— 

12.0 

0.0 

7 2 

4.0 

3.3 


3.3 x 10' 

I ,ovv 

— 

1 1 0 

S.7 

0 I 

3.S 

2.3 





Fraction ol Maximum 

lielativ e Lntrnpv 




Relat iv e 



Im Indicated 

Oils li.md 





Scene 







Sensor 


C .’omplevitv 

S 

7 

0 

3 

4 

3 

TM 


High 

l no 

0 03 

0.73 

0 33 

0.37 

0.27 



Medium 

1 no 

0. S3 

0 02 

0 1 1 

0 30 

0 21 



1 ,< >vv 

LOO 

0.77 

0 30 

0 3S 

0.21 

0.17 

MSS 


High 



1 .00 

0 70 

0.00 

0 40 

0.31 



Medium 

— 

L00 

0 S3 

0 00 

0 41 

0.28 



1 .ovv 

— 

1 no 

0.70 

0 33 

0.33 

0.22 
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tl ic* number of variables than oi tlu* t\ pc of variables 
(original bands or transformed). TM bad greater rel- 
ative entropy values Ibr Brightness and Brightness/ 
Greenness than did MSS. Information in the Tas- 
seled Cap Third Component of TM was mueli 
greater than that of MSS. both by itself and in com- 
bination with Brightness or Greenness, confirming 
TM s greater dimensionality. In ductions in the 
number ol bits used to encode data in each channel 
decreased the information content, affecting TM 
data proportionately more than MSS data so that, 
with live bits or less per band, the information in 
comparable sets was equal. 
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